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conditions.  Analysis  of  the  restoring  function  revealed  the  existence 
of  lower  and  upper  static  as  well  as  dynamic  burning  instability. 
Nonlinear  static  and  dynamic  stability  boundaries  were  determined  and 
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tube  are  discussed. 
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CHAPTER  1 

SCOPE,  OBJECTIVES  AND  BACKGROUND 

Sec.  1.1  - SCOPE  OF  THIS  WORK 

Sec.  1.2  - OBJECTIVES  AND  PLAN  OF  PRESENTATION 

Sec.  1.3  - HISTORICAL  BACKGROUND  ON  THE  THEORY  OF  NONLINEAR 

DYNAMIC  STABILITY  OF  SOLID  PROPELLANTS  COMBUSTION. 
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Sec.  1.1  - SCOPE  OF  THIS  WORK 

The  ultimate  purpose  of  this  work  is  to  formulate  the  ba- 
sic laws  governing  the  nonlinear  dynamics  of  solid  propellant 
burning  within  the  framework  of  a thermal  theory.  This  implies 
♦ that  : 

(1)  the  findings  should  be  readily  extended  to  other  cases 
of  heterogeneous  combustion  ; in  particular  the  results 
are  not  dependent  on  the  nature  of  solid  propellant. 

(2)  new  flame  models  are  not  being  proposed,  rather  a method 
is  offered  to  judge  how  good  any  unsteady  flame  model  is. 

(3)  information  regarding  steady  state  solutions  is  found 
as  a particular  case. 

(4)  interactions  combustion/t luidodynamics  (e.g.,  in  a ro- 
cket combustion  chamber)  are  not  considered. 

By  making  the  usual  set  of  assumptions  (see  Sec.  2.1),  one 
is  reduced  to  study  the  stability  properties  of  the  nonlinear 
heat  conduction  equation  in  the  condensed  phase  (the  standard 
Fourier  partial  differential  equation  with  time-varying  boundary 
conditions) . This  is  done  by  an  approximate  analytical  technique, 
but  keeping  the  nonlinearity  of  the  problem.  It  will  be  shown 
that  a nonlinear  algebraic  function,  called  static  restoring  fun- 
ction, can  be  determined  defining  nonlinear  static  and  dynamic 
burning  stability  boundaries  for  finite  size  disturbances.  The 
dynamic  stability  boundaries  are,  in  principle,  restricted  to 
forcing  functions  monotonically  changing  in  time  (a  common  case 
in  actuality) , but  are  valid  for  any  kind  of  transient  burning 
(e.g.,  depressurization  or  deradiation).  * 

It  should  be  explicetelyjkbserved  that  no  attempt  is  made  to 
explain  or  predict  the  steady  flame  structure  of  burning  propel- 
lants. This  is  supposed  to  be  assigned  in  terms  of  burning  rate 
vs  pressure  (experimental  data)  and  flame  temperature  vs  pressure 
(thermochemical  equilibrium  computations) . This  information,  in 
addition  to  a proper  flame  model,  is  enough  to  define  the  funda- 
mental nonlinear  burning  stability  properties  of  a given  propellant. 

It  is  emphasized  that  this  stability  analysis  can  be  applied 
to  any  flame  model  and  solid  propellant  composition,  but  the  ac- 
tual numerical  values  do  depend  on  the  specific  flame  model  and 
solid  propellant  composition  chosen.  The  results  can  be  verified 
by  computer  simulation  (i.e.,  by  direct  integration  of  the  basic 
set  of  governing  equations)  and/or  by  experiments.  While  plenty 
of  computer  results  are  already  available,  the  experimental  part 
of  this  work  is  entering  just  now  in  its  crucial  phase. 

Since  this  is  the  first  detailed  report  on  the  subject  after 
the  Ph.D.  thesis  of  one  of  the  authors  (see  below) , the  opportunity 
is  taken  to  summarize  most  of  the  findings  of  the  nonlinear  stabi-  » 


--  — - 


lity  analysis  under  examination  and  to  make  the  point  about 
the  status  of  the  theory.  To  avoid  a too  long  discussion,  the 
review  will  be  limited  to  the  common  case  of  static  restoring 
function  of  the  type  A-B-C  (see  Sec  2.5).  The  less  common, 
but  specially  interesting,  case  of  static  restoring  function 
of  the  type  A-B-C-D-E  will  be  reviewed  in  the  next  annual  re- 
port; only  some  hints  will  be  given  in  this  report.  The  com- 
plete literature  so  far  published  by  this  working  group  on  the 
topics  treated  in  this  report  is  listed  as  Refs.  1-10;  the  most 
exhaustive  reference  is  Ref.  3 which  was  recently  distributed. 
Reading  of  this  report,  complemented  with  some  portions  (later 
specified)  of  Ref.  3 and  with  Ref.  10,  should  give  a complete 
picture  of  the  situation. 
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Sec.  1.2  - OBJECTIVES  AND  PLAN  OF  PRESENTATION 

Consider  the  physical  system  of  Fig.  la  representing  a 
strand  of  solid  propellant,  subjected  to  a radiant  flux  impin- 
ging with  istantaneous ‘intensity  ( 1 — r ) ♦ I (t)  at  its  surface, 
burning  with  instantaneous  rate  R(t)  in  a closed  vessel  at 
instantaneous  pressure  P(t)  and  ambient  temperature  Tq.  The 
overall  problem  of  transient  solid  propellant  combustion  con- 
sists of  predicting  the  burning  rate  history  in  time. 

Specific  objectives  of  this  study  are  : 

(1)  distinction  between  static  and  dynamic  burning  regimes. 

(2)  definition  and  determination  of  nonlinear  static  stabi- 
lity boundaries. 

(3)  definition  and  determination  of  nonlinear  dynamic  sta- 
bility boundaries. 

(4)  prediction  of  dynamic  extinction. 

(5)  prediction  of  number  and  nature  of  static  solutions. 

(6)  influence  of  the  implemented  flame  model  on  the  above 
predictions . 

(7)  influence  of  the  integral  method  on  the  above  predictions. 

(8)  is  there  a self-sustained  oscillating  burning  regime? 

(9)  is  there  an  ignition  boundary? 

(10)  experimental  verification  of  the  above  predictions. 

Most  of  these  questions  will  be  dealt  with  in  what  follows. 
First,  a literature  survey  of  the  problem  is  offered  in  Sec. 1.3 
(this  literature  survey  is  taken  from  Ref.  3 and  is  reported  for 
the  convenience  of  the  interested  reader).  Then,  the  basic  set 
of  assumptions  and  equations  is  illustrated  in  Sec.  2.2  ; the 
transformation  of  the  PDE  into  an  approximately  equivalent  ODE  is 
performed  in  Sec.  2.4  ; properties  of  the  resulting  nonlinear 
algebraic  restoring  function  are  discussed  in  Sec.  2.5.  Static 
and  dynamic  burning  stability  are  respectively  treated  in  Sec. 

2.6  and  Sec.  2.7.  The  point  on  the  status  of  the  theory  is  made 
in  Sec.  2.8.  Experimental  problems  are  discussed  in  Sec.  2.9 
and  Ch.  3.  Conclusions  and  suggestions  for  future  work  are  col- 
lected in  Ch.  4. 
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Sec .1.3-  HISTORICAL  BACKGROUND  ON  THE  THEORY  OF  NONLINEAR 
DYNAMIC  STABILITY  OF  SOLID  PROPELLANT  COMBUSTION 

The  question  of  the  dynamic  extinction  by  depressuri- 
zation and  in  general  of  the  burning  stability  has  been  well 
debated  in  the  technical  literature  (Refs.  11-19),  but  few 
works  have  been  really  constructive.  The  erroneous  appli- 
cation of  the  quasi-steady  gas  phase  assumption  (Ref. 20), 
the  strong  limitations  due  to  the  use  of  linearized  theories 
and  the  empirical  nature  of  several  of  the  proposed  criteria 
(Refs.  12-13)  are  the  most  serious  drawbacks  in  this  area. 

For  a critical  review,  the  interested  reader  might  wish  to 
consult  the  work  by  Merkle  (Ref.  21b).  A paper  by  T'ien  (see 
below)  is  the  only  one  aiming  directly  at  establishing  an 
extinction  criterion  for  fast  depressurization.  Merkle  him- 
self did  not  formulate  any  criterion,  but  just  picked  up  an 
empirical  value  of  surface  temperature  ( 0g  =0.43)  below 
which  chemical  reactions  are  considered  tosbe  too  weak  in 
order  to  sustain  the  deflagration  wave.  Likewise,  no  crite- 
rion has  been  formulated  so  far  for  the  case  of  dynamic  ex- 
tinction by  fast  deradiation  reported  by  the  Princeton  group 
(Refs.  3,  22,  23b). 

T'ien  argues  that  "heat  losses  are  the  mechanism  for 
both  static  and  dynamic  extinctions  of  solid  propellants" 

(Ref.  24) ; this  view  is  not  fully  shared  in  this  instance 
(see  Sec.  2.7).  In  any  event  T'ien  concludes  that  "for  depres- 
surization transients,  if  the  instantaneous  burning  rate  drops 
below  the  unstable  burning  rate  solution  at  the  final  pressu- 
re, extinction  will  occur".  T'ien  derives  his  quantitative 
criterion  from  another  study  by  him  of  flammability  limits 
of  premixed  flames  under  the  influence  of  environmental  di- 
sturbances (Ref.  25) . A somewhat  similar  result  will  be  obtai- 
ned here  by  following  a completely  different  approach. 

The  line  of  research  evolved  within  the  framework  of 
the  mechanistic  (Zeldovich)  approach  has  been  unable  to  reach 
truly  meaningful  conclusions  so  far  about  the  dynamic  stabi- 
lity boundary.  Istratov  et  al.  in  1964  (Ref.  26)  used  an  in- 
tegral method  in  order  to  determine  an  approximate  solution 
to  the  unsteady  nonlinear  energy  equation  in  the  condensed 
phase  of  a propellant  burning  with  constant  surface  tempera- 
ture. Extinction  was  assumed  to  occur  when  the  surface  ther- 
mal gradient  on  the  condensed  side  exceeds  a critical  value 
corresponding  to  the  static  stability  line.  This  is  a serious 
mistake,  since  nothing  can  really  be  said  a priori  about  the 
possibility  of  a dynamic  burning  regime  in  a range  of  burning 
rate  that  is  statically  unstable.  In  this  study  also  an  inte- 
gral method  has  been  used,  but  completely  different  results 
have  been  found.  Novozhilov  in  1967  (Ref.  27)  improved  the 
previous  model  by  considering  a variable  surface  temperature 
and  recognizing  that  dynamic  burning  is  allowed  also  in  the 
range  of  parameters  where  statically  stable  solutions  are  not 
found.  Extinction  is  then  assumed  to  occur  when  the  burning 


r 
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rate  at  the  final  pressure  drops_  below  a limiting  value  ex- 
perimentally established  in  nonstationary  burning  conditions. 
This  criterion  heavily  relies  in  very  delicate  experimental 
results  and,  moreover,  fails  to  assign  any  premonitory  si- 
gnal, if  any,  about  the  approaching  of  the  limit. 

This  question  of  a possible  early  warning  of  extinction 
during  a depressurization  transient  evaluated  via  the  Zeldo- 
vich  method  was  examined  in  a paper  by  the  Princeton  group 

• (Ref.  28)  in  1972.  The  possibility  was  checked  that  the  cros- 
sing of  the  static  stability  boundary  is  sufficient  to  sub- 
sequently produce  dynamic  extinction.  No  clear  answer  to  this 
question  was  given.  The  position  taken  here  is  that  the  static 
stability  boundary  has  only  secondary  relevance  in  a dynamic 
situation.  Indeed,  according  to  the  same  Princeton  reference, 
"the  dynamic  conditions  of  extinguishement  tend  to  shift  the 

• stability  line"  (p.  257  of  Ref.  2,8)  . More  extensive  work  in 
this  area  (Ref.  29)  also  failed  to  reach  any  significant  con- 
clusion. In  any  event,  Novozhilov  (p.  216  of  Ref.  30)  in 
1973  observes  that  this  "question  requires  certain  informa- 
tion about  the  properties  of  the  system  outside  the  area  of 
smooth  burning.  Such  information  cannot  be  obtained  from 
experiments  on  steady  state  combustion.  For  the  calculation 
of  unsteady  conditions  in  the  unstable  region  it  is  necessa- 
ry to  draw  on  certain  schemes  of  combustion,  which  make  it 
possible  to  predict  the  properties  of  propellants  beyond 

the  (static)  stability  limit".  In  other  words  the  author  ad- 
mits the  failure  of  the  mechanistic  approach  in  this  respect 
and  recognizes  that  burning  in  the  statically  unstable  region 
can  only  be  treated  with  a flame  model.  However,  the  need  for 
more  advanced  treatments  in  the  framework  of  flame  models  is 
advocated  in  Ref.  31. 

In  this  work,  quantitative  criteria  for  dynamic  stabili- 
ty of  burning  propellants  (and  in  particular  for  dynamic  ex- 
tinction) are  defined  by  means  of  flame  models.  The  integral 
method  of  Goodman  (Ref.  32)  is  implemented  in  order  to  apply 
known  mathematical  methods  to  the  resulting  approximate  ODE 
formulation  of  the  problem.  The  same  method  has  already  been 
applied  (Ref.  33  )t  droplet  burning;  somewhat  simpler  use  has 
been  made  at  Princeton  (Refs.  34-35)  in  solid  propellant  ro- 
cket engines.  The  concept  of  using  the  simpler  ODE  formula- 
tion of  the  problem,  instead  of  the  PDE  one,  is  relatively 
common  in  Soviet  literature.  However,  the  method  of  transfor- 
mation is  rather  different;  an  interesting  review  of  the 
mathematical  problem  was  made  by  Gostintsev  (Ref.  36) . The 
approach  has  been  applied  mainly  to  ignition  (Refs.  37-38,  e.g.), 
unsteady  burning  (Refs.  26-39,  e.g.)  and  stability  problems 
(Refs.  30,40,  e.g.).  Remark  that  the  Soviet  approach  (Refs. 
26,30,36-40)  differs  not  only  in  the  mathematical  details 
but  especially  in  the  structure  of  the  physical  model.  Indeed, 
instead  of  using  a flame  model,  the  Soviet  investigators  re- 
sort to  the  Zeldovich  method,  consisting  of  constructing 


the  thermal  gradient  at  the  condensed  phase  side  of  the 
burning  surface  from  experimental  steady  state  data  (see 
Ref.  28,e.g.).  However,  it  was  already  mentioned  that 
the  Zeldovich  method,  requiring  steady  state  experiments , 
is  useful  in  establishing  intrinsic  stability  boundaries 
but,  at  least  in  principle,  cannot  be  extended  to  obtain 
dynamic  stability  boundaries. 
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In  this  chapter  an  overview  of  the  theoretical  deve- 
lopments, both  analytical  and  numerical,  concerning  nonli- 
near burning  stability  of  solid  propellants  is  offered. 
First  the  basic  assumptions  and  equations  are  revised, 
then  the  transformation  of  the  governing  PDE  into  an  ODE 
is  performed.  This  allows  us  to  define  static  and  dynamic 
stability  properties  of  a burning  solid  propellant.  For 
the  first  time: 

(a)  results  obtained  by  using  MTS  and  KTSS  flame  models 
are  compared; 

(fc)  the  effect  of  the  power  law,  chosen  for  implementing 
an  integral  method,  on  the  shape  of  the  static  resto- 
ring function  is  shown? 

(c)  the  nonlinear  static  stability  boundary  is  determined 
directly  from  the  flame  model  and  also  through  the  sta- 
tic restoring  function. 

Moreover,  some  detailed  analyses  of  the  MTS  and  KTSS  (both 
linearized  and  full  versions)  flame  models  are  offered.  Nu- 
merical computations  verifying  the  analytical  predictions 
are  shown.  A section  on  the  state  of  the  art  of  the  nonli- 
near burning  stability  analysis  ends  the  theoretical  deve- 
lopments of  this  chapter.  Experimental  results  from  the 
shock  tube  are  reported  in  the  last  section. 


Sec.  2.1  - BACKGROUND  AND  NOMENCLATURE 


The  physical  system  dealt  with  is  presented  in  its 
most  general  form  in  the  schematic  of  Fig.  la.  The  oressure 
of  the  vessel,  the  radiant  flux  impinging  on  the  surface  of 
the  strand  (originating  exclusively  from  some  external  ener- 
gy source) , the  ambient  temperature  measured  at  the  cold 
boundary  of  the  propellant  sample  (supposed  to  be  infinite- 
ly long) , and  any  other  parameter  which  can  be  controlled 
in  a known  way  from  the  outside  of  the  vessel  are  designa- 
ted as  controlling  parameters.  A change  of  one  or  more  of 
those  controlling  parameters  will  affect,  in  some  way,  the 
'•late  of  the  physical  system  and,  consequently,  they  are 
also  called  external  perturbations. 

On  the  other  side,  we  designate  as  intrinsic  perturba- 
tions sources  all  .those  "small"  (in  a sense  to  be  specified) 
irregularities  and  imperfections  always  present  in  the  real 
world. but  which  are  nevertheless  neglected  in  the  idealized 
picture  of  Fig.  la.  For  example,  nonuniform  composition 
of  the  strand,  impurities  variuosly  scattered  in  the  conden- 
sed phase,  complicated  geometry  of  the  burning  surface,  etc. 
all  contribute  to  hopefully  small  but  persistent  differences 
of  the  actual  phenomena  from  those  described  by  mathematical 
models . 

Whether  the  perturbation  sources  are  external  or  intrin 
sic,  the  prior,  supposedly  unperturbed  steady  state  profile 
of  temperature  in  the  combustion  wave  will  be  modified  to 
some  new  perturbed  unsteady  profile.  Let  us  define  the  di- 
sturbance temperature  profile  as  the  profile  of  the  point- 
by-point  difference  between  the  perturbed  profile  and  the 
original,  unperturbed  profile.  The  general  problem  of  static 
stability  may  be  stated  as  follows:  given  a stationary  sta- 
te of  the  physical  system,  the  system  is  forced  to  a close 
but  nonstationary  state  and  it  is  asked  whether  the  system, 
after  a long  enough  period  of  time,  will  go  back  to  its  ini- 
tial state  or  will  move  away  from  it.  In  the  specific  physi- 
cal configuration  considered' in  this  study,  it  is  asked  whe- 
ther the  disturbance  temperature  profile  will  die  out  in 
time  or  not.  Mathematically,  the  problem  is  an  initial  va- 
lue problem  and  is  usually  described  by  a parabolic  type  of 
partial  differential  equation. 

It  is  of  concern  to  distinguish  between  static  and  dy- 
namic stability.  The  general  problem  of  dynamic  stability 
may  be  stated  as  follows:  given  an  initial  stationary  state 
of  the  physical-  system,  the  system  is  forced  to  a different 
(final)  stationary  state  by  means  of  appropriate  changes  in 
time  of  pressure  and/or  radiant  flux  and  it  is  asked  whether 
the  system,  after  a long  enough  period  of  time,  will  reach 
the  wanted  final  state  or  another  (unwanted)  final  state.  In 
other  words,  the  stability  of  a system  where  intrinsic  per- 
turbation sources  exclusively  are  considered  to  be  acting 
(stability  of  a state)  is  called  static .Conversely , the  sta- 
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bility  of  a system  in  which  the  external  controlling  para- 
meters are  changing  in  time  (stability  of  a transition)  is 
called  dynamic.  It  will  be  shown  that  in.  the  latter  case, 
the  rate  of  change  in  time  of  the  pressure,  for  example, 
is  of  fundamental  importance  and  this  explains  the  expres- 
sion "dynamic  stability". 

In  this  research  program,  the  ultimate  objective  is  to 
establish  boundaries  separating  regions  of  stability  from 
those  of  instability  or.  some  convenient  graphical  plot  . It 
will  be  shown  that  static  and  dynamic  stability  boundaries 
are  different  in  nature. 

Since  the  mathematical  problem  is  formidable,  no  gene- 
ral method  has  been  found  so  far  for  solving  the  stability 
problem  in  its  various  aspects.  Historically,  a large  amount 
of  work  has  been  devoted  to  stability  problems  in  fluid  dy- 
namics. A standard  treatment  in  this  and  related  fields  is 
the  linearized  approach,  which  is  based  o<in  two  essential  as- 
sumptions: disturbance  quantities  infinitesimally  small  and 
a mathematical  model  (differential  equations,  boundary  con- 
ditions, and  any  other  relationship)  containing  only  distur- 
bance quantities  of  first  order.  Under  these  hypotheses,  the 
problem  reduces  to  an  eigenvalue  formulation  whose  mathema- 
tical theory  is  very  well  developed.  Besides  this  standard 
approach  , several  other  approximate  and  often  ad  hoc  methods 
have  been  set  up  in  various  fields  of  applied  sciences.  The 
reason  for  such  a confusing  state  of  affairs  is  essentially  < 

one:  the  behavior  of  most  physical  system  is  described  by 
nonlinear  equations  and  nowadays  not  only  a mathematical 
theory  encompassing  all  types  of  nonlineaarities , but  even  a- 
nalytical  methods  capable  of  dealing  with  specific  types  of 
nonlinearities,  are  not  yet  available. 

In  this  study,  in  order  to  avoid  the  serious  drawbacks 
of  any  linearized  treatment,  a nonlinear  combustion  stability 
analysis  by  means  of  an  approximate  approach  is  offered.  Ba- 
sically, the  problem  will  be  transformed  ifrom  a PDE  into  an 
ODE  formulation  via  an  integral  method.  Notice  that  the  pro- 
posed approach  requires  the  implementation  of  a flame  model, 
since  Zeldovich  method  is  of  no  help  (see  Sec.  1.3).  However, 
any  flame  model  capable  to  describe  quantitatively  the  heat 
feedback  from  the  gas  phase  to  the  burning?  surface  is  in 
principle  acceptable.  The  choice  of  the  OLame  model  is  a du- 
ty of  the  investigator.  Obviously,  the  proposed  stability 
theory  can  only  extract  the  stability  properties  (both  static 
and  dynamic)  implicitely  buried  in  any  flame  model.  This  fact, 
in  turn,  might  help  in  the  choice  of  the  fc>est  flame  model 
for  a given  combustion  problem. 

A composite,  ammonium  perchlorate  bashed  solid  rocket 
propellant,  denoted  as  AP/PBAA  No.  941,  is;  taken  as  datum  ca- 
se throughout  this  report.  However,  this  boy  no  means  is  meant 
to  imply  that  the  approach  is  restricted  to  some  specific 
class  of  solid  propellants.  The  physical  and  chemical  proper- 
ties of  the  propellant  AP/PBAA  No.  941  are  listed  in  Tab.1, 
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Sec.  2.2  - FORMULATION  OF  THE  GENERAL  QUASI-STEADY  GAS 
PHASE  TRANSIENT  PROBLEM 

Except  where  explicitely  excluded,  the  following  set 
of  assumptions  is  valid  throughout  this  work.  VJith  referen- 
ce to  Fig . 1 : 

a)  Entire  Strand. 

. Monodimensional. 

2.  Constant  thermal,  physical,  and  chemical  properties. 

3.  Condensed  phase  and  gas  phase  have  identical  specific 
heats . 

4.  At  cold  boundary,  in  thermal  equilibrium  with  ambient. 

5.  Subjected  to  no  external  forces  (acceleration,  gravity, 
electromagnetism) . 

6.  No  emission  of  radiation  (only  external  radiation  sour- 
ces) . 

b)  Condensed  Phase. 

1 . Semi-infinite  slab. 

2.  Uniform  and  isotropic  composition. 

3.  Adiabatic,  except  at  the  burning  surface. 

4.  No  chemical  activity. 

5.  No  radiation  scattering. 

6.  No  photochemical  effects. 

c)  Interface. 

1 . Infinitesimally  thin  plane  surface  (collapsed  reacting 
layer) . 

2.  One-way,  irreversible  gasification  process  (pyrolysis). 

d)  Gas  Phase. 

1.  Semi-infinite  column  of  gas. 

2.  Mixture  of  thermally  perfect  gases  of  average  molecu- 
lar weight  m. 

3.  One-phase,  laminar,  nonviscous,  strongly  subsonic  flow. 

4.  Adiabatic,  except  at  the  burning  surface. 

5.  No  interaction  with  (external)  radiation. 

6.  Lewis  number  = 1,  each  chemical  species  has  the  same 
specific  heat  C„.  mass  diffusion  is  expressed  by  Pick's 
law;  therefore  the  gas  phase  can  be  described  by  a sim- 
ple thermal  model  (see  Sec.  2.3). 

7.  Quasi-steady  behavior. 

The  following  set  of  reference  parameters  is  used  for 
nondimensionalizing  (with  specific  reference  to  the  propel- 
lant AP/PBAA  No.  941  taken  as  datum  case;  see  also  Tab.  1): 
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Pref 

= 

68 

atm 

Tref 

m 

300 

K 

Rref 

= 

R {pref } = °-837 

cm/s 

T 

s,ref 

Ts(Pref)=  1'000 

K 

Tf ,ref 

= 

Tf(Pref)=  2'430 

K 

xref 

5 

**C/Rref  = 1-673  x ID’3 

cm 

fcref 

= 

^c/Rref  = 1 *998  x 10"3 

s 

Qref 

= 

Cc(Ts,ref-Tref)  = 231 

cal/g 

Xref 

= 

?cCcRref(Ts,ref-Tref)  = 297  * 8 

cal/cm 

from  which  one  gets  the  following  nondimensional  variables: 


9 
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x 

<5 

a 

X 

H 

Q 

F_ 


= P/P 
= R/R 
= x/x 


ref 

ref 


ref 
(1/a)  /x 


ref 


= t/t 


ref 


“ VQref 
= Qf/Qref 
- Io/Iref 

“ i^Iref 


nondim.  pressure 
nondim.  burning  rate 
nondim.  distance 

nondim.  absorption  layer  thickness 
nondim.  time 

nondim.  surface  heat  release 
nondim.  gas  heat  release 
nondim.  radiant  flux  intensity 
nondim.  conductive  heat  flux. 


Notice  that  for  the  temperature,  both  of  the  following  defi- 
nitions are  used  according  to  need: 


(2.2.1a) 
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TW  - T~l 
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And  the  choice  of  the  particular  definition  depends  on  the 
specific  physical  phenomenon  being  considered. 

The  general  nondimensional  quasi-steady  gas  phase  tran- 
sient formulation  includes  the  unsteady  condensed  phase  pro- 
blem (Eq.  2.2.2),  relationships  accounting  for  the  surface 
mass  production  (Eqs.  2.2.3),  the  quasi-steady  gas  phase  pro- 
blem (Eqs.  2.2.4  - 2.2.8),  some  auxiliary  relationships  assi- 
gning the  gas  phase  and  burning  surface  heat  releases  (Eqs. 
2.2.9  - 2.2.10)  and  the  time  history  of  the  controlling 
parameters  (Eqs.  2.2.11  - 2.2.12).  Obviously,  the  gas  phase 
treatment  is  not  complete  without  the  energy  equation  whose 
first  integral  gives  the  heat  feedback  law;  this  is  discus- 
sed in  Sec.  2.3.  For  details  about  the  writing  of  the  equa- 
tions, the  reader  is  referred  to  Ref.  3. 
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where  the  boundary  condition  BC1  at  X=0  merely  states  the 
energy  conservation  at  the  burning  surface  (see  Fig.  1b). 

In  the  above  BC1 : 

is  the  nondimensional  heat  feedback 
to  be  defined  by  a flame  model; 

is  a nondimensional  heat  loss  from 
the  burning  surface; 

H is  conventionally  supposed  positive  if  endothermic; 

N is  a trasparency  factor  (N  =0  for  opaque  condensed  pha- 
se and  = 1 for  transparent  condensed  phase) . 

The  last  point  deserves  some  comment.  The  hypothesis  of 
a collapsed  reacting  layer  might  imply  in  some  cases  that  a 
certain  thickness  of  condensed  phase  Is  shrunken  to  an  infi- 
nitesimally small  value.  For  consistency,  then,  the  (external) 
radiant  flux  should  be  separated  in  a portion  deposited  direc- 
tly at  the  surface  (as  if  the  propellant  were  opaque)  and  in 
a (remaining)  portion  distributed  volumetrically  in  the  con- 
densed phase  according  to  the  proper  extinction  coefficient. 


The  trasparency  factor  Nfcf  therefore,  should  be  set  equal 
to  the  fraction  of  radiant  flux  absorbed  in  the  reacting 
layer  when  it  is  thought  of  as  extended.  However,  the  ex- 
tent and  the  very  existence  of  a condensed  phase  reacting 
layer  is  an  extremely  controversial  question.  In  this  stu- 
dy N is  taken  either  as  0 or  1 , although  this  is  not  a 
necessity. 


As  to  the  surface  mass  production,  the  following  com- 
bined pyrolysis  law  is  used  according  to  need  both  for  stea- 
dy and  unsteady  states: 


(2.2.3a) 


Arrhenius  pyrolysis,  &s 


(2.2.3b)  31  - 


KTSS  pyrolysis,  &H  » 0 > 


(2.2.3c) 


51  x 0 


where©  is  some  empirically  defined  minimum  temperature 
under  which  the  solid  propellant  undergoes  no  chemical 
activity.  The  power  w is  determined  by  matching  the  two 
pyrolysis  laws  at  a given  surface  temperature,  say©  = 0^, 
near  the  ambient  temperature.  This  formulation  of  tne 
pyrolysis  law  allows  us  to  avoid  the  "cold  boundary  diffi- 
culty" of  the  Arrhenius  expression. 

Relationships  for  the  gas  phase  include  the  mass  con- 
tinuity 


(2.2.4)  U(x,r)- ^(t)  ; 

(A  V 

the  momentum  equation  describing  a trivial  pressure  distri- 
bution. 


(2.2.5)  (X,  t)  - 


the  state  equation  (in  dimensional  terms) 
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the  integral  energy  balance  assigning  the  quasi-steady  flame 
temperature 


where  ^as  to  be  defined  by  a flame  model;  and  the 
integral  energy  balance  assigning  the  steady  state  flame  tem- 
perature in  presence  of  an  external  radiant  flux 


(2.2.8) 


. 0-^.-  it 

nfy.F.) 


where  (q  fCc  =1  has  been  assumed  and  O j is  the  adiaba- 

tic flame  temperature  furnished  by  standard  thermochemical 
computations.  For  example,  for  the  datum  case  considered 
(see  Tab.  1),  the  following  relationship  was  found  to  hold 
up  to  100  atm:  in  dimensional  terms 


Up) 


£ ,«£ 


(P*f  - p) 


50 
68 

= 68  atm  (Ref.  21b) . 

The  heat  release  in  the  gas  phase  is  considered  to  be 
uniquely  defined  by  the  pressure:  f or C . 1 one  finds 


where  T^  = 2430  K and  Pref 


(2.2.9) 


Q(9) 


w ♦ ©i  + 


©r 


The  heat  release  at  the  burning  surface  is  considered 
strictly  depend  on  the  original  condensed  phase  composition 
only;  for  example,  for  the  broad  class  of  composite,  AP  ba- 
sed propellants  fRef.  43c) 


. p (Qcr  * Qv,  AP  - Qa/pa  ) + (l-  Qvt  B , 

(2.2.10)  H « = 

Finally,  it  is  understood  that 

(2.2.11)  (y)  = externally  assigned 

(2.2.12)  ("t j = externally  assigned. 
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Sec.  2.3  - MTS  and  KTSS  FLAME  MODELS 

In  order  to  define  explicitely  the  function  s 50 
required  for  the  BC1  of  the  PDE  describing  the  trans-ient 
thermal  profile  in  the  condensed  phase  (Eq.  2.2.2),  a flame 
model  has  to  be  chosen.  It  is  emphasized  that  any  flame  mo- 
del can  be  picked  up.  For  the  composite  propellant  AP/PBAA 
No.  941  taken  as  a datum  case,  two  flame  models  were  tested: 
the  MTS  (Ref. 21)  and  the  KTSS  (Ref. 42)  flame  models.  Both 
were  developed  by  Summerfield  and  coworkers;  other  flame  mo- 
dels will  be  tested  in  the  near  future. 

The  model  worked  out  by  Merkle,  Turk,  and  Summerfield 
(MTS)  is  basically  a judicious  combination,  based  on  experi- 
mental data,  of  the  diffusive  and  kinetic  aspects  of  the  aas 
phase  reactions.  The  model  worked  out  by  Krier,  T'ien,  Siri- 
gnano,  and  Summerfield  (KTSS)  can  be  see  as  a particular  ca- 
se of  the  previous  one  when  the  flame  is  purely  diffusion 
controlled.  Both  models  are  thermal  in  nature  and  invoke 
the  assumptions  mentioned  in  Sec.  2.2.  The  MTS  model  will  be 
used  for  most  of  the  theoretical  developments  of  this  report. 
For  a detailed  knowledge  of  how  the  model  is  derived,  the 
reader  is  referred  to  the  original  Ref.  21b  and  also  to  Refs. 
43-44,  treating  the  GDF  theory  from  which  the  MTS  flame  model 
has  evolved.  Only  qualitative  comments  will  be  offered  here. 

Any  pattern  of  chemical  reactions,  no  matter  how  compli- 
cated, involves  the  two  physically  distinct  phenomena  of  mass 
diffusion  and  chemical  kinetics.  It  is  known  from  homogeneous 
combustion  theory  that  quite  often  diffusion  processes  domi- 
nate the  kinetic  processes  in  the  sense  that 


T 


di 


» 1 


or,  conversely,  the  opposite  extreme  may  be  true.  These  two 
limiting  configurations  are  respectively  known  as  diffusion 
and  premixed  flames.  For  heterogeneous  combustion,  in  parti- 
cular of  solid  propellants,  there  are  instances  in  which 
the  gas  phase  portions  of  the  deflagration  wave  may  be  trea- 
ted according  to  one  of  the  above  limiting  configurations. 
But,  in  principle,  especially  for  transients  connected  to 
large  excursions  of  the  controlling  parameters,  the  gas  pha- 
se treatment  of  a heterogeneous  combustion  wave  has  to  ac- 
count for  both  processes.  This  was  originally  done  by  Sum- 
merfield in  his  GDF  theory  for  the  AP  steady  flame  model.  La- 
ter, an  extension  was  made  for  the  AP  quasi-steady  flame  mo- 
del accounting  for  variable  flame  and  surface  temperatures. 
Although  derived  for  the  specific  case  of  the  AP  class  of 
propellants,  the  formulation  of  the  MTS  flame  model  is  quite 
general  and  may  be  adapted  to  other  cases  of  heterogeneous 
combustion. 


For  the  specific  case  of  an  AP  composite  propellants, 
the  MTS  flame  model  accounts  not  only  for  the  granular  dif- 
fusion flame  (GDF)  of  the  oxidizer  with  the  binder  (Ref.  21) 
but  also  for  a distended  premixed  flame  corresponding  to  the 
AP  decomposition  products  burning.  In  this  study,  the  simpli- 
fied version  with  the  premixed  flame  collapsed  at  the  bur- 
ning surface  has  been  adopted.  This  assumption  is  acceptable 
(p.  28  of  Ref.  21b)  for  pressure  values  not  lower  than  say, 

1 atm.  The  most  important  parameter  characterizing  the  gas 
phase  during  a quasi-steady  transient  is  the  heat  feedback 
to  the  burning  surface.  According  to  the  (collapsed)  MTS 
flame  model,  the  nondimensional  heat  feedback  is 


(2.3.1)  ^ = Q ft 


where  1^.  is  a nondimensional  kinetic  time  parameter; 

is  a nondimensional  diffusion  time  parameter; 

t1  is  a nondimensional  reaction  time  parameter; 
re 

% is  given  by  an  Arrehenius-type  pyrolis  law. 
Following  the  MTS  flame  model  development,  we  put 
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where  a second  order  gas  phase  reaction  has  been  postulated 
to  occur  wholly  at  the  highest  temperature  Tf  (p.29  of  Ref. 
43a)  and 


where  the  diffusional  mixing  rate  of  fuel  pockets  with  the 
surrounding  atmosphere  of  oxygen-rich  gases  is  assumed  to 
depend  on  the  surface  temperature  T (p . 3 1 of  Ref.  43a). 

S 

The  two  constants  A„  and  B„  are  determined  for  each 
specific  propellant  by  "the  best  fit  of  the  steady  state 
burning  rate  theory  to  the  measured  burning  rate  data"  (p. 
38  of  Ref.  21b).  An  application  is  shown  in  Fig.  2 for  the 
propellant  AP/P3AA  No.  941 . In  the  spirit  of  the  GDF  theo- 
ry, it  is  further  assumed  (e.g.,  see  p.  33  of  Ref.  21b) 
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that  the  overall  solid  propellant  reaction  time  parameter 
can  be  expressed  as  the  following  simple  combination  of 
the  two  above  limiting  cases: 


(2.3.4) 


This  relationship  has  been  shown  to  represent  the  pressure 
dependence  of  the  burning  rate  quite  accurately  for  a wide 
range  of  composite  solid  propellants  (Ref. 43)  and  also  to 
describe  the  depressurization  extinction  correctly  for  both 
composite  and  double-base  solid  propellants  (Ref. 21).  Howe- 
ver, it  should  be  recognized  that  the  whole  MTS  approach 
depends  on  the  arbitrary  assumption  of  Eq.  2.3.4  and  the 
choice  of  the  constants.  How  this  choice  is  affected  by 
the  fitting  procedure  mentioned  above  is  illustrated  in 
Fig.  2 (different  pressure  interval)  and  Fig.  3 (different 
low  end  of  the  pressure  interval) . 

According  to  the  KTSS  flame  model,  the  nondimensional 
heat  feedback  is 


(2.3.5.) 


% $ 


where  x'  is  an  average  nondimensional  diffusion  time  para 
meter; 

31  is  given  by  a power  pyrolysis  law. 

Notice  that  the  above  KTSS  heat  feedback  law,  derived  for  a 
dif fusionally  controlled  AP  flame,  is  a particular  case  of 


the  MTS  heat 

feedback 

law 

where 

(2.3.6) 

>>jX  K- 

(2.3.7) 
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That  is,  when  the  diffusion  time  parameter  is  much  lar- 
ger than  the  kinetic  time  parameter  and  is  temperatu  - 
re  independent.  Extensive  experimental  investigation 
( Ref.  43  ) shows  that  for  AP  propellant  flames  the  dif- 
fusional  mixing  is,  indeed,  sensibly  slower  than  the 


J 


chemical  kinetics  as  long  as  the  pressure  is  larger  than 
say,  5 atm.  It  is  also  expected  (Ref. 21)  that  flame  tem- 
perature undershoots  during  decelerated  transients  slow 
down  the  chemical  kinetics  much  more  effectively  than  the 
dif fusional  mixing  (compare  Eq.  2.3.2  vs  Eq.  2.3.3). 

In  its  linearized  version,  the  nondimensional  KTSS 
heat  feedback  law  is  given  by 


(2.3.8)  Q 

^'S  % Tj. 

which  is  a particular  case  of  the  full  expression  (Eq.2.3.5), 
when  the  further  assumption  is  made  that 

(2.3.9)  e « ± 


corresponding  to  tj;  very  large  but  finite  burning  rate. The 
constraint  of  Eq.  2.3.9  is  usually  valid  for  steady  burning, 
but  is  not  acceptable  in  extinction  transients,  since  under 
this  circumstance  'S.  — O while  remains  finite  and,  therefore, 


e 


1.0 


Note: in  the  linearized  version  commonly  used,  the  nondimen- 
sional KTSS  heat  feedback  law  is  written  as 


(2.3.9) 


?2B  ( 9 w) 


Then , bv comparing  Eq.  2.3.8  with  Eq.  2.3.9,  the  average  non- 
dimensional  diffusion  time  parameter  can  be  evaluated 


cp 

(2.3.10)  * ..  * 

9"  w)  •/ 

An  essential  ingredient  for  both  MTS  and  KTSS  flame  mo- 
dels is  the  experimental  stationary  31  (?)curve.  The  steady 
state  structure  of  the  MTS  and  KTSS  models  are  respectively 
shown  in  Fig.  4 and  Fig.  5 for  the  propellant  AP/PBAA  No. 941. 
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It  is  worthwhile  to  remark  that  even  the  linearized  version 
of  the  KTSS_model  can  accurately  reproduce  the  experimental 
stationary  ft  curve . A comparison  of  the  two  plots  shows 
no  strong  difference  as  to  the  steady  state  behavior  of  the 
flame  models  under  examination.  Consider  now  the  quasi-steady 
behavior  as  plotted  in  Fig.  6 (MTS  model)  and  Fig.  7 (KTSS 
model).  It  is  obvious  that  the  linearized  KTSS  model  is  phy- 
sically unrealistic  for  low  burning  rate. 

In  summary: 

(a)  stationary  or  small  perturbation  solutions  of  the  MTS 
and  KISS  (even  linearized)  flame  model  are  similar  at 
high  pressure,  since  in  these  instances  only  small  chan- 
ges in  temperature  and  burning  rate  are  considered; 

(b)  at  low  pressure  and/or  low  burning  rate,  AP  flame  is  ki- 
netically  controlled  and  therefore  should  not  be  descri- 
bed by  the  KTSS  (even  non  linearized)  model; 

(c)  the  KTSS  linearized  model  can  never  predict  extinction 
kscause  the  heat  feedback  law  is  found  to  be  inversely 
proportional  to  the  burning  rate. 


A useful  comparison  for  what  follows  is  shown  in  Fig.  8, 
where  the  MTS  heat  feedback  (often  used  in  this  study)  and 
the  linearized  KTSS  heat  feedback  (commonly  used  in  the  li- 
terature) are  plotted  together  in  the  range  10  to  40  atm 
( <^0*  curves).  From  this  graphical  plot  one  observes  that,  for 
a given  set  of  external  conditions,  two  steady  solutions 
(reacting  configuration  A and  unreacting  configuration  C) 
are  found  for  the  MTS  model,  but  only  one  for  the  linearized 
KTSS  model  (reacting  configuration  A) . 


The  quasi -steady  gas  phase  working  map  for  the  MTS  flame 
is  illustrated  in  Fig.  9.  On  this  plot,  besides  the  usual  heat 
feedback  lines  at  constant  pressure,  lines  at  constant  ambient 
temperature  (for  adiabatic  burning)  or  at  constant  radiant  flux 
(at  standard  ambient  temperature)  are  drawn  measuring  the  sta- 
tionary heat  flux  


;*)-  $(©,-  8„)  * ^H-  F. 


absorbed  in  the  condensed  phase.  Where  crosses 

j ^ ) * steady  state  solutions  “ft  ( *9t  ^ A|«* 

are  singled  out.  Obviously,  the  curves  ^ n°t 

depend  on  the  choice  of  the  flame  model,  while  the  curves 
q.  , ('?,%)  do  depend.  Similar  remarks  hold  true  in  Fig.  10, 
whi're  a large  fan  of  lines  at  constant  ambient  temperabre  are 
drawn  on  a MTS  flame  plot.  Notice  that  even  at  zero  degrees 
Kelvin  (assuming  the  flame  model  still  valid  ! ) the  steady 
solution  is  far  from  the  lower  static  stability  boundary  (cf. 
Sec.  2.6). 


It  will  be  shown  that  the  burning  stability  properties 

are  strongly  affected  by  the  value  of  the  surface  heat  release. 

For  our  datum  case  (see  Tab.  1),  the  best  value  computed  from 

the  literature  is  Q = -158.2  cal/g  (the  exothermicity  is 

s 
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due  to  the  primary  ammonia/perchloric  acid  flame  occurring 
in  the  ammonium  perchlorate  decomposition) . Changing  Q , 
while  Keeping  ^(P|  fixed,  means  to  consider  a familys 
cf,soli<3  propellants  whose  R (P)  behavior  is  shown  in  Fig. 
11  and  Tab.  2. 
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Sec.  2.4  - AN  ORDINARY  DIFFERENTIAL  EQUATION  FORMULATION 
OF  THE  PROBLEM 

Basically,  the  mathematical  method  is  the  one  set  up 
for  the  first  time  by  Von  Karman  and  Polhausen  in  the  stu- 
dy of  boundary  layers  and  later  generalized  by  Goodman, 
among  others,  to  a large  number  of  thermal  problems.  The 
method  can  be  extended  to  any  other  problem  described  by 
nonlinear  PDE  of  parabolic  type.  In  our  case,  the  approach 
consist  s in  defining  a parameter  ) ( z ) , called  the  pene- 
tration distance  of  the  thermal  wave  in  the  condensed  pha- 
se, "such  that  for  (X)  the  propellant  slab,  for  all 

practical  purposes,  is  at  an  equilibrium  temperature  and 
there  is  no  heat  transferred  beyond  this  point"  (p.53  of 
Kef.  32) . The  evolution  in  time  of  the  thermal  profile  in 
the  condensed  phase  is  obtained  by  following  the  time  hi- 
story of  the  penetration  distance  propagating  into  an  ini- 
tially uniform  temperature  field.  Within  this  penetration 
layer,  progressing  in  time,  the  qualitative  space  distribu- 
tion of  temperature  is  assumed  known  a priori;  but  we  ma- 
ke sure  that,  in  so  doing,  the  integral  balance  of  thermal 
energy  in  the  condensed  phase  is  preserved . In  other  words, 
the  price  to  be  paid  for  making  the  transformation  from 
PDE  to  ODE  formulation  is  an  approximate  solution  of  the 
local  space  distribution  of  temperature;  this  is  not  such 
a serious  drawback  because  the  interest  is  in  the  time  evolu- 
tion of  the  surface  temperature.  In  any  event,  several  inve- 
stigations found  an  error  of  only  some  percents  for  various 
cases  in  which  both  the  exact  and  the  integral  solutions 
were  evaluated.  For  example,  see  Fig.  5 on  p.  89  of  Ref.  32 
showing  the  temperature  time  history  at  the  surface  of  a 
semi-infinite  slab  with  triangular  surface  heat  flux. 

In  order  to  get  a deeper  understanding  of  questions 
related  to  the  integral  method  as  applied  to  thermal  pro- 
blems, the  interested  reader  might  wish  to  consult  the  ex- 
cellent review  by  Goodman  (Ref. 32)  and  the  references  given 
there.  Before  getting  involved  in  mathematical  details,  the 
reader  should  be  warned  about  the  limits  of  the  integral 
method.  Any  solution  obtained  by  the  integral  method  con- 
tains hopefully  small  but  irrevocable  errors  in  the  final 
numerical  results.  The  question  of  how  to  improve  the  accu- 
racy then  arises.  It  has  been  argued  that  "there  is  no  a 
priori  guarantee  that  increasing  the  order  of  the  polynomial 
(used  to  represent  the  space  distribution  of  the  unsteady 
temperature  profile)  will  improve  the  accuracy.  Although 
the  accuracy  is  frequently  improved  with  this  technique, 
it  can  be  demonstrated,  nonetheless,  that  there  are  cases 
for  which  it  actually  worsens"  (p.96  of  Ref.  32).  In  this 
same  reference  it  is  suggested  that  the  method  of  weighted 
residuals  provides  a very  efficient  method  for  improving 
the  accuracy  of  the  results  obtained  by  using  the  integral 
method. 
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The  simple  integral  method  implemented  in  this  work 
can  be  considered  a special  case  of  the  method  of  weighted 
residuals  when  just  one  parameter  (the  penetration  distance 
£)and  one  weighting  function  are  considered.  Under  these 
circumstances,  the  only  equation  to  be  considered  is  the 
heat  balance  integral  (see  below) . Since  the  integral  method 
assumes  a uniform  initial  distribution  of  temperature,  a n ew 
nondimensional  variable  is  defined: 


(2.4.1) 


where 


(2.4.2) 


* (x,T)  • et(xj-  e(x,t) 


gi  (x)  - (©,  - ej  . 


X%L 


is  the  initial  steady  state  distribution  of  temperature 
and  9(X*) is  the  temperature  distribution  following  some 
perturbation.  The  new  variable  u(X,t)  may  therefore  be  conve- 
niently interpreted  as  the  finite  temperature  disturbance 
propagating  inside  the  condensed  phase  and  superimposed 
on  the  initial  temperature  distribution  after  the  action  of 
perturbation.  At  the  initial  instant  X =0,  by  definition 
e(3(,T)=  9;(x)and  u^r)=o,  Suppose  now  that,  in  the  following 
instant,  a perturbation  starts  acting  on  the  system  and  ma- 
kes u(X,r)^0:  the  goal  of  the  analysis  is  to  determine  the 
ultimate  effect  of  such  a temperature  disturbance  after 
waiting  a period  of  time  sufficiently  long  for  the  perturba- 
tion to  disappear.  No  assumption  whatsoever  is  made  as  to 
the  size  of  the  temperature  disturbance. 

The  analysis  will  be  restricted  to  the  case  of  an  op- 
tically opaque  propellant  (in  the  sense  that&/G)/(l/a)>>  i)at 
constant  (not  necessarily  reference) ambient  temperature 
The  basic  set  of  nondimensional  equations  for  an  optically 
opaque  propellant  (N  =0) , initially  burning  with  constant 
rate  31; at  constant  pressure  , while  subjected  to  a radiant 
flux  of  constant  intensity  ; , is 

f CIO  /o  0©  I Pt>£ 
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where  $ (x)  is  the  burning  rate  defined  by  any  appropriate 
pyrolysis  law  and  ) (QO/OX’)  ^ A is  the  heat  feedback 

law  defined  by  any  appropriate  flame  model.  In  terms  of 
the  disturbance  temperature  r),  the  set  of  equations  is: 


(2.4.4) 


where 


(2.4.5) 


and  ( A-^/Ac)  (3u/3X)q  5 is  again  defined  by  any  appropriate 
flame  model.  ' 

Notice  that  at  the  initial  instant  T=0,  the  surface  energy 
balance  is 


(2.4.6) 


where 


49; 


\ l 


K\AXj  3, 

E. 


s -H 


(2.4.7)  < 


31; . e! 


1 

/ ^r«f 

1 - 

k T-F 

\kF 

- i 


for  ©:  s * ©K 


for  .<©* 


TBlS  TkQt  ^ jo  DPfl 


r 


26 


Therefore, 


(2.4.8) 


The  transformation  into  an  ODE  of  the  PDE  of  Eq.2.4.4, 
with  BC1  given  by  Eq.  2.4.8,  can  be  performed  by  assuming 
a polynomial  dependence  of  the  disturbance  temperature  on 
the  space  variable: 

(2.4.9)  u'CX'r')  “ K°  Cx)*  \ (r)X  * K3  ft)X  - - 


BC1  of  Eq.  2.4.4  can  be  written  as 


Comments  about  the  implications  of  this  particular  thermal 
profile  are  given  below.  The  above  n+1  coefficients  k;  (r) 
are  to  be  determined  from  the  boundary  conditions  which  ex- 
press no  disturbance  (up  to  the  (n-2)-th  derivative)  at  the 
cold  end  of  the  penetration  depth  and  the  energy  balance  at 
the  hot  boundary  of  the  penetration  depth.  With  our  formula- 
tion this  implies: 


(2.4.10) 
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where  is  the  penetration  distance  to  be  obtained  be- 

low. After  algebric  manipulations,  one  obtains  for  the  unknown 
profile  of  the  disturbance  temperature 


For  X = 0 one  obtains  at  the  burning  surface 
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from  which  the  unknown  penetration  depth  can  be  expressed 


(2.4.11) 


The  disturbance  temperature  profile  is  then  given  by 


(2.4.12)  u(X/'rj  « 


where  u (T)  is  the  unknown  surface  temperature  disturbance 
to  be  determined.  It  is  obvious  from  the  above  relationship 
that  the  time  history  of  the  disturbance  temperature  is  re- 
stricted to  disturbance  thermal  profiles  monotonically  de- 
caying in  space . Indeed , there  is  no  way  for  a polynomial 
profile  to  accommodate  an  inflection  point.  This  difficulty 
may  be  overcome  by  "allowing  a new  penetration  depth  to  be- 
gin propagating  at  each  maximum  or  minimum"  (p.96  of  Ref. 32) 
of  the  heat  input  into  the  condensed  phase.  However,  this 
procedure  will  not  prove  necessary  here,  since  testing  both 
the  static  and  dynamic  stability  does  not  require  a detailed 
knowledge  of  the  structure  of  the  thermal  profile. 

A space  integration  can  now  be  performed  over  X from 
0 to  - l(^)  of  the  PDE : 


(2.4.13) 


where 


(2.4.14) 


-JM 


being  defined  by  Eq.  2.4.11.  Uoon  substituting  the  approxi- 
mate disturbance  temperature  profile  described  by  Eq.  2.4.12,  the 
following  differential  expression  is  obtained: 
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which  is  approximately  equivalent  to  the  initial  PDE+boundary 
conditions  formulation.  Notice  that  in  the  above  equation 
©;  £X) is  known  (see  Eq.  2.4.2)  and  that  only  the  surface  tem- 
perature appears. 


The  above  nonlinear  ODE  in  the  unknown  u£  (T)  describes  the 
response  of  the  system  to  a finite  size  departure  of  the 
surface  temperature  from  the  stationary  value  due  not  only 
to  intrinsic  perturbation  sources  acting  on  the  system  (sta- 
tic stability) , but  also  to  any  arbitrary  but  externally 
assigned  monotonical  change  in  time  of  controlling  parameters 
such  as  pressure  and  radiant  flux  (dynamic  stability) . 

If  one  wishes  to  know  the  v/hole  temoerature  profile,  he 
has  only  to  substitute  u (t)  in  Eq.  2.4.12,  evaluate!  (r) 
from  Eq.  2.4.11  with  the  help  of  Eq.  2.4.8,  and  then  use  the 
u MO  definition  (Eq.  2.4.1)  in  order  to  determinate  the 
resultant  temperature  profile. 


Considerations  of  a general  character  on  the  static  sta- 
bility of  the  system  described  by  Eq.  2.4.19  can  now  be  made. 
According  to  the  theory  based  on  the  Lyapunov  theorems,  a 
given  equilibrium  configuration  of  the  system  is  asymptotical- 
ly stable  if 

S(T)-**o  1 


j 

In  other  words,  the  system  is  (asymptotically)  stable  is  the 
disturbance  disappears  at  large  time  and  the  system  returns 
to  its  initial  configuration.  Following  Lyapunov  (p.  216  of 
Ref.  45),  Eq.  2.4.19  can  be  written,  in  a more  concise  way,  as 
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dependent) 

If  no  forcing  function  is  acting  on  the  system,  i.e. 
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cost  and 


then  Eq.  2.4.20 


reduces  to 


F = cost, 
o 


(2.4.23) 


by  which  the  autonomous  contribution  assumes  the  meaning 
of  a static  restoring  function.  Indeed,  under  these  circum- 
stances, one  can  think  of  the  chemically  reacting  system 
comprised  of  the  deflagrating  condensed  material  as  analo- 
gous to  a mass-spring  type  of  mechanical  system  with  the 
nonlinear  characteristic  f (u  ) . Considering  the  u s (t)  defi- 
nition of  Eq.  2.4.1,  Eq.  2.4:20  can  also  be  written  as 

• / 


The  static  restoring  function  f ~ depends  only  on 

the  nature  of  the  deflagrating  substance;  the  nonautonomous 
term  changes  in  time  according  to  the  first  derivatives  of 
the  external  controlling  parameters  (see  Eq.  2.4.22). 


Sec.  2.5  - NATURE  OF  THE  STATIC  RESTORING  FUNCTION 

The  static  restoring  function  is  an  algebraic  nonlinear 
function  strictly  dependent  on  the  nature  of  the  burning  pro- 
pellant (including  the  specific  way  its  flame  is  described,  i.e. 
the  flame  model  implemented  and  the  order  of  the  approxima- 
ting polynomial  used  in  Sec.  2.4)  and  the  operating  conditions. 
Before  examining  quantitative  plots,  consider  the  qualitative 
plots  of  Fig.  12a  and  Fig.  12b.  This  allows  to  extract  the  basic 
properties  of  the  static  restoring  function  independenly  on  the 
flame  description.  It  is  anticipated  that  the  heat  release  at 
the  surface  of  the  propellant  is  the  most  infuential  parameter 
on  the  shape  of  the  static  restoring  function. 

Consider  the  qualitative  plot  of  Fig. 12a.  According  to 
Eq.  2.4.24,  when  no  forcing  function  is  acting_on  the  system, 
all  the  points  (algebraic  roots)  for  which  f (9t/i  - ©4  ) =0 

define  possible  equilibrium  configurations  for  the  burning 
propellant,  since  they  correspond  to  =0.  It  is  seen 

in  Fig.12athat,  in  addition  to  the  trivial  ©s  =0  (unburning 
propellant,  root  C) , two  more  equilibrium  solutions  (roots 
A and  B)  are  allowed,  in  general,  for  a given  set  of  parame- 
ters. Let  us  consider  the  equilibrium  configuration  corre- 
sponding to  root  A.  Suppose  that,  for  some  unspecified  rea- 
son, the  burning  rate  or  the  surface  temperature  increases 
a finite  amount;  then  the  burning  prooe_llant  is  no  longer 
in  a stationary  configuration , dQ^jdx  «=  - ©4)  is  nega- 

tive, and  the  system  reacts  by  decreasing  its  surface  tempe- 
rature. Conversely,  if  for  an  unspecified  reason,  the  surface 
temperature  of  a propellant  burning  according  to  the  confi- 
guration of  root  A decreases,  the  system  reacts  by  increa- 
sing its  surface  temperature.  These  movements,  around  point 
A,  are  indicated  by  arrows  in  Fig. 12a.  It  is  concluded  that 
the  equilibrium  configuration  corresponding  to  root  A is 
stable  because,  when  disturbed,  the  system  always  moves  back 
toward  A.  By  the  same  arguments,  it  is  concluded  that  the 
equilibrium  configuration  corresponding  to  root  B is  unsta- 
ble; any  disturbance  yields  movement  away  from  the  point. 
Therefore,  in  steady  state  experiments,  only  solution  A is 
observed . 

If  the  set  of  parameters  is  changed,  for  example  the  pres- 
sure is  increased  from  P to  P. , the  initial  condition  of  the 
ODE  (see  Eq.  2.4.2)  is  changed  so  that  a new  pair  of  roots,  A. 
and  B^ , is  found  in  addition  to  the  trivial  ©s  = 0 of  root  C. 
Again,  root  A^  defines  a stable  equilibrium  configuration, 
while  root  defines  an  unstable  equilibrium  configuration 
corresponding  to  the  new  set  of  parameters.  Likewise,  a new 
pair  of  roots,  A2  and  B_,  in  addition  to  the  trivial  ©5  =0  of 
root  C,  is  found  if  the  pressure  is  decreased  from  P to  P~. 

It  follows  that  the  axis  in  Fig .1 2a  includes , in  addition 
to  the  trivial  solution  of  nonburning  propellant  at  the  root 
C,  a segment  of  stable  solution  A.  and  a segment  of  associa- 
ted unstable  solution  B.  (each  pair  of  roots  corresponds  to 
a given  set  of  parameters) . 
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It  should  be  explicitely  observed  that  the  trivial  =0 
solution  nay  be  obtained  only  if  the  exponential  Arrhenius-type 
pyrolysis  law  for  mass  production  at  the  surface  is  dropped. 
Obviously,  the  trivial  0S  = 0 solution  implies  that  no  external 
energy  source  (e.g.,  radiation)  is  acting  on  the  solid  propellant. 

The  qualitative  picture  of  Fig.  12a  illustrates  the  general 
behavior  of  the  static  restoring  function  when  the  pressure  is 
varied  parametrically  at  fixed  ambient  temperature  and  surface 
heat  release.  The  behavior  of  the  static  restoring  function  when 
the  surface  heat  release  is  varied  parametrically,  at  fixed  am- 
bient temperature  and  pressure,  is  illustrated  by  the  qualitati- 
ve picture  of  Fig. 12b.  Of  course,  all  previous  considerations 
hold  true.  For  example,  when  the  surface  heat  release  is  low 
enough  (in  a sense  which  will  be  better  understood  below) , the 
system  behaves  again  according  to  the  static  restoring  function 
represented  by  curve  CBA  (Fig.  12b).  However,  for  increasing  va-  ' 
lues  of  the  surface  heat  release,  it  is  found  that  the  static 
restoring  function  is  represented  by  curve  CB.A.D.E. . This  is 
rather  surprising;  in  principle,  there  are  now  five  equilibrium 
configurations . 

Based  on  the  previous  analysis,  however,  it  is  immediately 
recognized  that  C is  the  stable  equilibrium  solution  for  the  un- 
reacting state  (trivial  solution) . But  further  discussion  is 
required  to  understand  the  nature  of  the  remaining  four  roots. 

This  is  easily  accomplished  by  considering  the  steady  state  e- 
nergy  balance  of  the  overall  combustion  wave.  For  example,  the 
graphic  plot  of  Fig.  9 suggests  that  for  each  set  of  parameters 
only  one  solution  exists  for  the  reacting  state.  Conventionally, 
let  us  call  A that  particular  root  of  Figs.  12a  and  12b  corre- 
sponding to  the  energy  balance  solution  of  Fig.  9.  It  follows  that 
root  , although  stable  according  to  Lyapunov,  is  eliminated 
as  being  a reacting  equilibrium  solution.  In  other  words,  root 
E.  is  a false  equilibrium  solution  introduced  by  the  approximate 
ODE  formulation  of  the  problem.  The  remaing  roots  B^  and  D1  are, 
then,  both  unstable  equilibrium  solutions  for  the  reacting  sta- 
te. 

For  further  increase  of  the  surface  heat  release,  it  is  ob- 
served that  A-  and  D-type  roots  respectively  increase  (moving 
to  right)  and  decrease  (moving  to  left)  in  the  plot  of  Fig.  12b, 
until  coalescence  and  then  crossing  over  occur  with  exchange  of 
stability  character.  This  important  point  will  be  discussed  in 
detail  later  (see  Sec.  2.6).  For  further  increasing  of  the  surfa- 
ce heat  release,  B-  and  D-type  roots  disappear  after  coalescence, 
while  both  A_  and  E^  roots,  for  different  reasons,  are  eliminated 
as  being  stable  reacting  solutions  (curve  CA-E-  in  Fig.  12b). 

Under  these  circumstances,  it  follows  that  the  only  allowed  solu- 
tion is  the  trivial  unreacting  configurations  represented  by  root 
C.  Any  attempt  to  produce  a stationary  combustion  wave  with  a 
static  restoring  function  of  type  CA_E«  will  inevitably  result 
in  extinction.  This  type  of  extinction,  however,  cannot  be  quali- 
fied as  "dynamic". 
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_ Quantitative  plots  of  the  static  restoring  function 
$ (6n  - 0j)  vs  the  nondimensional  surface  temperature  0S  are  gi- 
ven for  the  propellant  AP/PEAA  Mo.  941.  A quantitative  plot 
requires  the  choice  of  a specific  flame  model  and  a specific 
order  of  the  approximating  polynomial  (see  Sec.  2.4).  By  ap- 
plying different  flame  models  to  the  same  propellant,  diffe- 
rent stability  properties  are  predicted  : this  offers  a cri- 
terium  for  discriminating  potentially  good  from  bad  flame  mo- 
dels. In  this  v:ork  MTS  and  KTSS  (both  linearized  and  non  li- 
• nearized)  flame  models  are  implemented.  As  to  the  order  of 

the  approximating  polynomial,  a cubic  lav/  v/as  chosen  to  repre- 
sent the  space  distribution  of  the  disturbance  thermal  profile. 
This  choice  was  suggested  by  a large  body  of  literature  on  heat 
transfer  problems  (e.g.,  Ref.  32)  and  by  similar  solid  propel- 
lant rocket  problems  (Refs.  33-34).  However,  there  is  no  a 
priori  guarantee  that  is  the  best  (and  even,  simply,  a good) 
choice.  The  choice  has  to  be  verified,  somehow.  This  is  shown 
in  Sec.  2.7. 

Plots  obtained  implementing  the  MTS  and  KTSS  flame  models, 
with  a cubic  disturbance  thermal  profile,  are  respectively  shown 
in  Figs.  13-16  and  in  Figs.  17-19.  The  MTS  flame  model  was  ap- 
plied by  assuming  a combined  exponential  power  pyrolysis  law 
at  the  burning  surface  (see  Eqs.  2.2.3,  with  0^=0. 15  and  0m  = O)  . 
The  KTSS  flame  model  was  applied  by  assuming  a power  pyrolysis 
law  (see  Eq.  2.2.3b,  with  w=6  and  =0)  over  the  whole  range 
of  surface  temperature  of  interest.  Comparative  results  are 
shown'  in  Fig.  20  for  some  standard  operating  conditions, 
leading  to  the  conclusion  that  MTS  in  perhaps  the  best  flame 
model  but  KTSS  nonlinearized  is  strongly  competitive.  The  ef- 
fect of  the  order  of  the  approximating  polynomial  for  the  di- 
sturbance thermal  profile  is  shown  in  Figs.  21-23  for  the  MTS 
flame  model  only. 

The  effect  of  the  ambient  pressure  on  the  static  restoring 
funption  vs  the  nondimensional  surface  temperature  is  illustra- 
ted in  Fig.  13  for  the  indicated  set  of  parameters  (MTS  flame 
and  n=3) . An  increase  of  pressure  implies  an  increase  of  the 
stable  equilibrium  surface  temperature,  but  a less  pronounced 
increase  of  the  unstable  equilibrium  surface  temperature.  The 
strength  of  the  stability  (see  Sec.  2.6)  is  enhanced  by  an  in- 
crease of  pressure.  Similar  comments  hold  true  as  to  the  effect 
of  an  external  radiant  flux  on  the  restoring  function  (Fig. 14). 
It  is  important,  however,  to  note  that  an  increase  of  the  ex- 
* ternal  radiant  flux  decreases  the  unstable  equilibrium  surface 

temperature  and  above  a certain  value  of  radiant  flux  intensity 
(e.g.,  40  cal/cm  -s  for  the  set  of  parameters  in  Fig.  14)  no 
more  unstable  solutions  are  found.  This  implies  that,  in  princi- 
ple (see  Sec.  2.7),  at  each  pressure  a minimum  value  of  exter- 
nal radiant  flux  intensity  exists  above  which  the  dynamic  boun- 
dary can  no  longer  be  defined.  The  effects  of  surface  heat  relea 
se  (Figs.  15a  and  15b)  are  more  involved.  Note  that  for  increa- 
sing values  of  the  surface  heat  release,  the  behavior  of  the 
static  restoring  function  shifts  from  a CBA  - type  (see  Fig. 15a) 
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curve  for  Q =-150  and  -158.2  cal/g,  to  CBADE  (see  Fig. 15b) 
for  Q = -1^0  and  -180  cal/g,  and  to  CBDAE  for  Q.  = -200  cal/g, 
with  S-D  roots  coalescence  occurring  at  Q =-190  cal/g.  This 
implies  that  at  30  atm  the  system  will  show  dynamic  stability 
effects  not  only  in  the  low  range  of  surface  temperature  (root 
B in  Fig.  15a)  but  also  in  the  large  range  (root  D in  Fig.  15b). 

The  effects  of  ambient  temperature  (Fig.  16)  are  of  little  in- 
terest in  the  range  investigated. 

Similar  trends  are  observed  for  the  nonlinearized  KTSS  fla- 
me model,  again  with  n=3.  The  effect  of  the  ambient  pressure 
in  Fig.  17  is  very  close  to  that  found  with  MTS  (cf.  with  Fig. 

13)  ; minor  differences  are  found  in  the  region  near  the  sta- 
tically unstable  root.  The  same  comments  hold  true  as  to  the 
effect  of  surface  heat  release  (Fig.  18). 

Drastically  different  results  are  observed  for  the  lineari- 
zed KTSS  flame  model,  alv/ays  with  n=3.  The  effect  of  the  ambient 
pressure  in  Fig.  1 9 is  virtually  the  same  (as  compared  both  to 
MTS  and  KTSS  nonlinearized)  for  surface  temperature  about  or  lar- 
ger than  the  steady  state  value,  but  is  physically  meaningless 
for  surface  temperaure  less  than  about  90%  of  the  steady  state 
value  (no  zero  solution  and  no  unstable  root) . 

The  static  restoring  functions  for  the  three  flames  are 
plotted  simultaneously  in  Fig.  20,  always  with  n=3 , for  the 
standard  conditions  specified  in  the  figure  (P=30  atm,  Ta=300  K, 

Q = -158.2  cal/g,  adiabatic).  The  comparison  graphically  em- 
phasizes the  differences  just  mentioned.  The  linearized  KTSS 
is  discarded  from  being  in  principle  an  acceptable  flame  model, 
the  nonlinearized  KTSS  is  acceptable  for  most  of  the  surface 
temperature  range  of  values  (except  near  extinction, since  che- 
mical kinetics  is  not  considered) , the  MTS  is  in  principle  accep- 
table over  the  whole  range  of  surface  temperature  of  interest. 

The  effect  of  the  order  of  the  polynomial  approximating 
the  disturbance  thermal  profile  on  the  static  restoring  function 
is  illustrated  in  Fig.  21  (P=20  atm),  Fig. 22  (P=30  atm)  and  Fig. 

23  (P=40  atm)  for  the  MTS  flame.  The  shape  of  the  static  restoring 
function  is  kept  for  n varying  from  2 to  4 ; the  steady  state  (or 
stable  reacting)  solution  is  always  recovered  (by  construction) ; 
the  trivial  zero  solution  is  also  alv/ays  recovered;  the  surface 
temperature  associated  with  the  unstable  reacting  solution  slightly 
increases  for  n decreasing.  Numerical  values  can  be  better  appre- 
ciated in  Tab.  3:  typically, the  surface  temperature  associated 
with  the  B root  increases  of  less  than  10%  when  n decreases  from 
3 to  2. 
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Sec . 2.6-  NONLINEAR  STATIC  STABILITY  ANALYSIS 

Static  stability  analysis  of  a burning  propellant  relates 
to  the  capability  of  the  propellant  to  keep  its  burning  equili- 
brium configuration  in  time.  It  was  shown  in  Sec.  2.4  that  this 
implies  to  study  the  properties  of  the  nonlinear  ODE 


which  in  turn  depends  on  the  shape  of  the  algebraic  nonlinear 
svatic  restoring  function  f (0,;s-  Os)  examined  in  Sec.  2.5.  The 
following  problems  are  of  relevance  in  the  area  of  static  sta- 
bility: 

1 . prediction  of  number  and  nature  of  the  allowed  static  solu- 
tions ( see  § 2.6.1). 

2.  prediction  of  static  stability  boundary  (see  § 2.6.2). 

3.  measurement  of  static  stability  strength  (see  £ 2.6.3). 

4.  prediction  of  the  pressure  deflagration  limit. 

The  first  three  problems  only  will  be  treated  in  this  report. 


§ 2.6.1  - Number  and  nature  of  the  allowed  static  solutions 

It  was  shown  in  Ref.  3 that  three  regions  of  interest  can 

be  distinguished,  according  to  the  value  of  Q ( see  Tab.  4) : 

s 

1.  before  A - D roots  coalescence,  the  static  reacting  solutions 
is  the  usual  stationary  combustion  wave  strictly  defined  by 
root  A. 

2.  after  A - D roots  coalescence  but  before  B - D roots  calescen- 
ce,  the  static  reacting  solution  is  a self-sustained  oscillating 
combustion  wave  around  root  A. 

3.  after  B - D roots  coalescence,  no  static  reacting  solutions 
whatsoever  is  allowed. 

Obviously,  the  trivial  unreacting  solution  0 is  al- 

ways allowed  (being  statically  stable).  Therefore,  too  much  sur- 
face energy  release  excludes  any  steady  reacting  solution  and 
under  these  circumstance^  the  only  possible  configuration,  even 
in  a static  environment, is  the  unreacting  state.  These  predictions 
were  verified  by  numerical  integration  of  the  governing  PDE. 

See  Figs.  24-25.  The  presence  of  self-sustained  oscillations  in 
Fig.  24  suqqests  the  existence  of  a limit  cycle.  Being  a limit 
cycle  an  overall  property  of  the  governing  partial  differential 
equation,  once  triggered  this  oscillatory  behavior  would  not  de- 
pend on  the  initial  conditions  of  the  system.  Indeed,  it  was 
found  in  Ref.  2 that,  independently  on  its  past  history,  under 
the  appropriate  circumstances  the  combustion  wave  undergoes  sharp 
self-sustained  oscillations  around  root  A with  amplitude  nearly 
defined  by  roots  D (lower  peak)  and  E (upper  peak)  and  with  pe- 
riod corresponding  to  the  characteristic  time  of  the  condensed  • 
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phase  thermal  wave. 

Under  the  standard  set  of  parameters  indicated,  numerical 
values  for  the  allowed  static  solutions  are  given  as  an  example 
in  Tab.  4.  In  particular,  the  range  of  Q values  for  v/hich  self- 
sus'tained  oscillating  combustion  waves  afe  expected  is  given, 
in  function  of  pressure,  in  Tab.  5.  Both  calculations  were  per- 
formed for  a MTS  flame  with  n=3. 

These  questions  will  further  be  investigated. 


$ 2.6.2  - The  static  stability  boundary 


An  important  piece  of  information  in  the  static  stability 
analysis  of  a burning  propellant  is  the  minimum  value  of  surfa- 
ce temperature  (or  burning  rate)  above  v/hich  stable  steady  bur- 
ning is  allowed.  In  order  to  evaluate  this  value,  a most  important 
feature  common  to  all  static  restoring  functions  is  emphasized. 

The  specific  surface  temperature  value,  0^,  at  which  the  pair  of 
solutions  Af  stable  and  Bf  unstable  coalesces,  at  a given  pressu- 
re (Fig.  12a)  and  for  a given  set  of  parameters, 


t 
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defines  the  branching  of  metastable  point  at  that  pressure  (and 
for  that  set  of  parameters) . Corresponding  to  this  special  value 
of  the  surface  temperature  a change  of  the  character  of  stabili- 
ty occurs  at  the  given  pressure,  in  that  for  ©5  < Gf  all  the  equi- 
librium solutions  at  that  pressure  are  statically  unstable  (roots 
B.)  while  for  ©5  > Of  all  the  equilibrium  solutions  at  that  pressu- 
re are  statically  stable  (roots  A.).  Therefore,  the  branching 
point  at  a given  pressure  isolate!  the  critical  static  stabili- 
ty value,  0/,of  the  surface  temperature  at  that  pressure.  Con- 
struction of  the  static  stability  boundary  on  the  *31  vs  ^pla- 
ne, in  a range  of  pressure,  consists  of  connecting  the  critical 
static  stability  points  defined  for  each  pressure  of  interest  in 
the  wanted  range  (and  for  a given  set  of  parameters).  This  implies 
the  search  of  the  branching  point  at  any  fixed  pressure. 


A general  method  (based  on  nonadiabaticity  of  the  propellant) 
for  solving  this  problem  is  discussed  below.  However,  for  the 
specific  case  of  a solid  composite  propellant , the  search  for  the 
branching  point  may  be  performed  immediately,  although  not  rigoro- 
usly,based  on  observation  of  the  physical  processes.  Consider  the 
qualitative  picture  superimposed  in  Fig.  26  on  a standard  MTS  flame 
model  plot.  Recall  that  burning  rate  and  surface  temperature  are 
univocally  related  by  an  appropriate  pyrolysis  law  . Therefore,  it 
is  irrelevant  to  define  the  branching  point  in  terms  of  burning 
rate  or  surface  temperature.  Now,  as  schematically  indicated  in 
Fig.  26,  the  branching  point  is  determined  at  each  pressure  by 
the  geometrical  condition  of  vertical  tangent,  .i.e. , infinite  slo- 
pe or  maximum  hesat  feedback.  Indeed,  for  0S  > ©f  a change  of  bur- 
ning rate  or  surface  temperature  produces  opposing  effects  on  the 
thermal  gradients  at  the  surface  from  the  condensed  and  the  gas 
phase  sides;  for  example,  the  former  increases  but  the  latter  de- 
creases when  the  burning  rate  is  disturbed  upward.  This  competi- 
tion is  the*  stabilizing  mechanism.  For  © < 0*  a change  of  burning 
rate  affects,  in  the  same  direction,  both  thermal  gradients;  this 
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is  destabilizing .Just  at  the  point  with  vertical  tangent,  the  two 
unstable  and  stable  static  roots  A = B coalesce  and  an  exchange 
of  the  stability  character  of  the  iolutfons  occurs;  therefore,  we  , 
may  identify  this  point  as  the  wanted  branching  or  metastable  point 
0 j .For  Qs  - Qf the  behavior  of  the  system  depends  on  the  direction 
of  the  perturbation.  Now,  in  order  to  construct  a static  stability 
boundary  in  the  plane  51  vs  q.  it  is  enough  to  connect  the  branching 
points  of  t 4.!*  sassociatcd  with ’each  pressure.  This  boundary  separates 
the  region  where  static  solutions  are  allowed  from  the  region  where 
no  static  solution  can  be  found. 

This  method,  first  suggested  in  Ref.  3,  is  not  rigorous  since 
the  contribution  from  the  surface  reactions  is  neglected.  Its  ad- 
vantage is  simplicity.  For  accurate  results,  the  procedure  illustra- 
ted in  Figs.  27  (Qo  - -120  cal/g) , Fig.  28  ( Q =-  158.2  cal/g)  and 
Fig.  29  (Q  =-200  cal/g)  for  a MTS  flame  has  tosbe  used.  This  consists 
in  assigning  a larger  and  larger  heat  loss  to  the  burning  propellant 
until  just  one  solution  is  found  (coalescing  from  the  stable  and 
unstable  static  roots) . Obviously,  this  search  can  be  made  directly 
for  any  flame  model  without  knowledge  of  the  static  restoring  func- 
tion. However,  the  same  results  should  be  obtained  by  searching  the 
coalescence  of  the  corresponding  A and  B roots  of  the  static  resto- 
ring function.  This  is  illustrated  in  Fig.  30  (P=20  atm),  Fig.  31 
(P=30  atm)  and  Fig.  32  (P=40  atm)  for  the  MTS  flame  with  n=3.  The 
closeness  of  the  two  set  of  results  confirms  the  validity  of  the 
ODE  formulation  of  the  problem.  Obviously,  this  second  approach  is 
affected  by  the  choice  of  the  approximating  polynomial  order.  The 
quantitative  differences  between  the  two  set  of  results  may  be  appre- 
ciated by  comparison  of  Tab.  G (heat  feedback)  with  Tab.  7 (static 
restoring  function  with  n=3)  . 

• 

£ 2.6.3  - Measurement  of  static  stability  .strength 

The  static  stability  analysis  can  be  put  on  a quantitative  ba-  t 
sis  by  means  of  the  "first"  Lyapunov  criterion.  According  to  this, 
the  solution  of  the  nonlinear  autonomous  ODE  of  Eq.  2.6.1  is  stable 
in  the  neighborhood  of  a point  ©*  if 

df  (^i,s  " ^s)  < O 

at  the  point  of  interest  0S  . The  physical  meaning  of  this  criterion 
was  explained  in  Sec.  2.5.  In  keeping  with  the  spirit  of  this  physi- 
cal interpretation,  the  strength  of  the  stability  of  the  solution 
at  a point  9*  may  be  measured  by  the  magnitude  of  the  derivative 
df/d  O evaluated  at  9,=  . How  this  value  is  affected  by  typical 

parameters  is  shown  in  Tab.  8.  In  the  exceptional  case  of  a bran- 
ching, or  metastable,  point  in  which  the  unstable  and  the  stable 
solutions  coalesce,  the  above  first  criterion  is  not  valid.  In  the 
present  situation,  however,  it  is  enough  to  say  that  the  point  ©f 
is  stable  from  one  side  (Qs^Qf+Z)  while  it  is  unstable  on  the  other, 
side  (©j  * Of  - e) , so  that,  overall,  the  point  9/is  said  to  be  unsta- 
ble. It  should  be  clear  that  the  first  Lyapunov  criterion  (used  for 
measuring  the  strength  of  the  stability)  is  essentially  a lineariza- 


tion  criterion  of  the  originally  nonlinear  ODE.  However,  it  is 
plajkible  that  both  roots  C and  A are  stable  to  finite  disturbances 
of  an  unknown  size,  also.  The  extent  of  the  stability  region  around 
the  stable  roots  could  be  estimated  by  the  Lyapunov  "direct"  method; 
but  this  is  not  one  of  the  purposes  of  this  analysis.  Physically, 
it  is  expected  that  the  range  of  stability  of  C root  (always  stable) 
is  limited  upward  by  B, while  the  range  of  stability  of  A root  (if 
stable)  is  limited  downward  by  B.  But  dynamic  effetcs  may  increase 
the  range  of  stability  of  C root  against  the  range  of  stability  of 
A root. 

Values  of  the  static  stability  strength  in  various  conditions 
are  summarized  in  Tab.  8. 


Sec.  2.7  - NONLINEAR  DYNAMIC  STABILITY  ANALYSIS 


It  is  wished  to  predict , under  an  appropriate  set  of 
assumptions,  the  behavior  of  an  heterogeneous  deflagration 
wave  initially  propagating  steadily  under  a given  set  of  ex- 
ternal parameters  and  then  subjected  to  arbitrary  but  knov/n 
changes  in  time  of  pressure  and/or  radiant  flux.  This  is  ea- 
sily done  if  an  ODE  formulation  of  the  problem  is  available. 

Dynamic  stability  of  a burning  propellant  relates  to 
transition  of  the  propellant  from  a given  initial  configura- 
tion to  a wanted  final  configuration.  A burning  propellant 
'is  dynamically  stable  if  the  transition  is  successful;  if  the 
transition  leads  to  a final  configuration  different  from  the 
wanted  one,  the  burning  propellant  is  dynamically  unstable. 

It  was  shown  in  Ref.  3 that  for  the  wide  range  of  con- 
trolling parameters  varying  in  time  according  to  monotonically 
decreasing  functions  (see  Fig.  33),  the  statically  unstable 
root  of  the  static  restoring  function  evaluated  at  the  final 
operating  conditions  is  also  the  asymptotic  (in  time)  dynamic 
limit  (see  Fig.  34).  It  was  verified  in  Ref. 4 that  this  limit 
holds  true  for  any  monotonic  decrease  in  time  of  pressure  or 
radiant  flux  intensity.  It  was  shown  in  Ref.  3 that  this  limit 
is  a unique  property  of  the  nropellant  at  the  final  operating 
conditions;  in  par ticular itdoes  not depend  on  optical  trasparency 
of  the  condensed  phase. 

Dynamic  effects  are  associated  both  with  a B-typo  root 
(see  Figs.  12a  - 12b)  and  with  a D-type  root  (see  Fig.  12b). 

The  former  is  called  lower  dynamic  stability  and  relates  to 
extinction;  the  latter  is  called  upper  dynamic  stability  and 
relates  to  vigorous  accelerations  of  the  cor.bustion  wave  pos- 
sibly followed  by  extinction  due  to  cverstabil  ty  (Fig.  24) . 
Discussions  and  results  concerning  this  poin  are  given  in 
Refs.  2-3. 

Note  that  in  the  general  ODE  formulation  of  Eq.  2.4.19 
the  influence  of  time  varying  external  parameters  is  felt 
only  through  their  time  derivatives.  This  implies  that,  if  no 
change  of  pressure  and/or  radiant  flux  occurs,  the  system  is 
only  subjected  to  random  intrinsic  fluctuations  and  the  results 
from  the  static  stability  analysis  apply.  It  is  also  implied 
that,  in  the  case  that  g (T , 9 . - ©o)  « f (6.  - 9^)  in 

Eq.  2. 4. 24,  the  dynamic  stabilit^sboun£lary  is  a1h§-return 
boundary  for  any  t(even  finite)  and  for  any  external  law  (even 
non  monotonic  or  increasing) . 

A graphical  summary  of  this  nonlinear  s Lability  analysis 
is  offered  in  Fig.  35  (MTS  flame  plot)  and  F g.  36  (standard  (R 
vs  P plot) . According  to  the  value  of  the  su  face  heat  release, 
upper  static  and  dynamic  stability  boundaries  may  also  be  found. 

Checks  of  the  lower  dynamic  stability  boundary  were  per- 
formed for  several  deradiation  and  depressurization  transients. 
See  Tabs.  9-10  and  Figs.  37  - 47.  The  computer  checks  were 
performed  by  solving  numerically  the  PDE  formulation  of  the 
general  quasi-steady  gas  phase  transient  problem,  whereas  the 
stability  boundaries  were  determined  by  analytical  means  from 


. 


considerations  of  the  approximately  equivalent  ODE  formulation. 
Since  the  lower  dynamic  stability  boundary  corresponds  to 
statically  unstable  roots,  it  cannot  be  observed  directly  either 
» experimentally  or  theoretically.  A go/no-go  technique  is  requi- 

red (see  Fig.  37).  The  results  shown  in  Figs.  38  - 47  and  summa- 
rized in  Tab.  10  suggest  that  n = 3 gives  an  excellent  agreement 
for  final  pressures  not  larger  than,  say  30  atm.  For  large  values 
of  pressure,  2 <n<3  seemes  more  successful  (see  Fig.  45). 


42 


Sec.  2.8  - STATE  OF  THE  ft RT  OF  THE  NONLINEAR  APPROACH 

The  problem  of  the  quasi-steady  qas  phase  burning  of  a so- 
lid propellant  in  a rather  general  form  is  being  dealt  with. 

In  order  to  retain  the  nonlinearity  of  the  problem,  an  aoproxi- 
mate  formulation  in  terns  of  an  ODE  was  written.  This  was  done 
by  means  of  an  integral  approach  limited  to  situations  in  which 
no  inflection  point  in  the  history  of  the  condensed  phase  thermal 
profiles  would  occur.  It  is  shown  that  the  important  facts  (see 
below)  of  heterogeneous  combustion  may  be  ascertained.  The  two 
fundamental  regimes  of  the  static  (intrinsic  random  perturbations) 
and  of  the  dynamic  (externally  assigned  changes  of  the  control- 
ling parameters)  stability  of  heterogeneous  flames ,allowing  for 
finite  size  disturbances , are  examined. 

The  following  facts  emerge  from  the  nonlinear  static  stabi- 
lity analysis.  For  a given  set  of  parameters: 

1 . a stable  stationary  nonreacting  equilibrium  configura- 
tion (trivial  solution)  is  always  found. 

2.  a stable,  stationary  (low  Qs)  or  self-sustained  oscilla- 
ting (large  Q ),  reacting  equilibrium  configuration  may 
be  found. 

3.  a stable  stationary  reacting  equilibrium  configuration 
is  found  before  A-D  roots  coalescence  occurs. 

4.  self-sustained  oscillations  are  found  between  A-D  roots 
coalescence  and  B-D  roots  coalescence. 

5.  after  B-D  roots  coalescence,  no  stable  reacting  equili- 
brium solution  is  found. 

6.  the  static  stability  boundary,  for  small  but  finite  si- 
ze random  intrinsic  disturbances,  is  the  locus  of  the 
branching  points  (A-B  roots  coalescence) . 

7.  the  stability  of  the  steady  state  equilibrium  configura- 
tions can  be  measured  by  the  slope  of  the  static  resto- 
ring function  vs  the  surface  temperature  at  the  point 

f5.  „ of  interest, 
l , s 

8.  The  effect  of  the  relevant  parameters  can  be  easily  eva- 
luated by  considering  the  corresponding  static  resto- 
ring function. 

The  following  facts  emerge  from  the  nonlinear  dynamic  sta- 
bility analysis.  For  a given  mcnotonic  law  of  timewise  decrease 
of  the  controlling  parameters: 

1 . extinction  nay  occur  even  though  the  final  point  of  the 
transition  is  statically  stable. 

2.  the  lower  dynamic  stability  boundary,  for  finite  size 
disturbances  consequent  to  timewise  monotonical  changes 
of  the  controlling  parameters,  is  the  locus  of  the  sta- 
tically unstable  roots  (B-type)  associated  with  the  fi- 
nal operating  conditions. 
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3.  “tflfe-  lower  dynamic  stability  boundary  holds  true  both 

for  deradiatl'on  and  depressurization  runs, and  is  not 
affected  by  the  specific  forcing  law. 

4.  upper  dynamic  instability  is  related  to  D root  and 
may  cause  vigorous  acceleration  of  the  combustion 
wave  or  dynamic  extinction. 

The  following  facts  emerge  by  comparing  several  flame  mo- 
dels: 

1 . the  KTSS  linearized  flame  model  is  physically  meanin- 
gless for  burning  rate  less  than  about  90%  of  the  stea- 
dy value. 

2.  the  KTSS  nonlinearized  flame  model  is  physically  que- 
stionable for  burning  rate  near  zero. 

3.  the  MTS  flame  model  is  in  principle  acceptable  over 
the  whole  range  of  burning  rate. 

The  following  facts  emerge  by  changing  the  order  n of  the 
polynomial  approximating  the  disturbance  thermal  profile: 

1.  C and  A roots  are  not  affected. 

2.  B,  D,  and  E roots  are  affected. 

3.  the  static  stability  boundary  can  be  evaluated  indepen- 
dently of  n. 

4.  the  dynamic  stability  boundaries  (both  lower  and  upper) 
are  determined  in  function  of  n. 

5.  the  polynomial  of  order  n=3  gives  accurate  results  up 
to  30  atm  of  final  pressure,  while  3>  n > 2 is  more  sui- 
table for  larger  values  of  pressure. 

The  physical  meaning  of  the  stability  boundaries  is  the 
following.  The  static  stability  boundary  defines  a line  of  se- 
paration between  a region  where  stable  steady  state  solutions 
are  allowed  and  a region  where  only  unstable  or  no  steady  solu- 
tions are  found.  The  dynamic  stability  boundary  defines  that 
ultimate  burning  condition  (in  terms  of  burning  rate  or  surface 
temperature)  beyond  which  extinction  necessarily  follows.  The 
fundamental  importance  of  this  distinction  is  stressed  by  the 
fact  that,  under  dynamic  conditions,  the  propellant  may  momen- 
tarily burn  also  in  a region  of  statically  unstable  configura- 
tions. It  is  shown  that  the  dynamic  stability  boundary  collap- 
ses to  the  static  stability  boundary  when  the  rate  of  change 
of  the  externally  controlled  parameters  is  negligible. 

Notice  that  numerical  values  have  been  given  only  for  a 
particular  composite  propellant  (AP/PBAA  No.  941).  This  was 
done  simply  because  the  properties  of  that  particular  propel- 
lant and  good  flame  models  were  readily  available.  It  is  felt, 
however,  that  all  analyses  were  conducted  from  a broad  point 
of  view,  and  i:i  no  way  were  they  dependent  on  the  particular 
type  of  propel Lant  chosen  as  the  datum  case.  Therefore,  the 
conceptual  results  are  expected  to  hold,  although  in  different 
ranges  of  the  relevant  parameters,  for  any  kind  of  solid  pro- 
pellant. 
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Sec.  2.9  - SHOCK  TUBE  EXPERIMENTS 


The  pressure  deflagration  limits  of  the  three  propellant 
formulations  under  consideration  (an  AP-based  composite,  a cata- 
lyzed DB  and  a noncatalyzed  DB)  were  determined  in  the  actual 
operating  conditions  of  the  shock  tube  (Fig.  48a) . These  expe- 
riments were  conducted  statically,  i.e.  a propellant  sample  was 
placed  in  the  test  section  of  the  shock  tube  (end  wall)  and  igni- 
tion was  attempted  by  means  of  an  electrically  heated  wire  (Fig.  48b) 
The  ignition  was  considered  successful  only  if  a self-sustained 
flame  would  develop.  The  test  was  repeated  for  several  pressurizing 
gases  (nitrogen, air , oxygen)  and  for  ambient  pressure  parametrical- 
ly changed.  Results  for  the  composite  propellant  are  shown  in 
Fig.  49  (not  self-sustained  flame  means  that  a visible  flame 
v.'.ild  develop  but  disappear  immediately  after  an  electrical  hca- 
’ iny  of  several  seconds  duration) . The  samples  were  cylindrical 
pellets  of  8 mm  diameter  and  about  1 0 mm  thickness.  Tests  with 
noncatalyzed  DB  samples  of  1 6 mm  diameter  did  not  show  any  appre- 
ciable difference. 

Samples  were  then  placed  at  the  same  location,  ignited  by  an 
electrically  heated  wire  and,  once  steadily  burning,  subjected  to 
shock  waves  of  different  strength.  The  light  emitted  by  the  bur- 
ning propellant  was  detected  by  means  of  a photodiode  directio- 
nally selective  and  sensitive  in  the  visible  and  near  infrared 
(spectral  response  peak  at  0.8  /tm)  . Several  tests  with 
different  propellants,  ignition  pressures,  pressurizing  gases,  and 
shock  speeds  failed  to  show  conclusive  trends  as  to  the  dynamic 
burning  rate  behavior  (as  revealed  by  the  light  emission) . The 
basic  difficulty  is  the  fast  succession  of  shock  waves  and  expan- 
sion fans  impinging  on  the  surface  of  the  burning  propellant  du- 
ring the  same  test.  For  example,  see  Figs.  50-51  (in  which  the 
period  of  the  pressure  pulse  is  about  10  ms).  In  this  series  of 

tests,  complete  extinguishment  of  the  propellant  was  never  achieved. 

• 

In  order  to  overcome  the  problem  just  mentioned,  the  shock 
tube  was  slightly  modified  to  work  as  a piston  tube.  At  present, 
a teflon  piston  weighting  150  g is  being  used.  In  this  operating 
configuration,  a longer  period  (about  50  ms)  of  the  pressure  pul- 
ses can  be  obtained  and,  in  general,  a better  control  of  the  whole 
experimental  apparatus  is  possible.  Work  in  this  direction  is  now 
in  progress.  Preliminary  results  obtained  with  a noncatalyzed  DB  are 
shown  in  Figs.  52-55.  Again  the  propellant  was  located  at  the  end 
walhof  the  shock  tube  and  electrically  ignited  before  the  run.  The 
test  of  Figs.  52-53  v/as  run  with  an  initial  pressure  of  0.5  atm 
of  air;  the  emitted  light  reached  the  maximum  value  at  the  peak 


pressure  and  then  decreased  during  the  expansion.  In  this  test 
the  propellant  was  extinguished  and  then  re  ignited  by  the  subse- 
quent shocks.  A test  conducted  at  an  initial  pressure  of  0.75 
atm  of  air  is  shown  in  Figs.  54-55.  In  this  case  the  emitted  light 
reached  the  maximum  value  before  the  peak  pressure  and  then  decre- 
ased while  the  pressure  was  still  rising.  In  this  test  the  propellent 
was  completely  extinguished. 
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NOMENCLATURE 


-1 


M 


B, 


M 


= volume  absorption  coefficient,  cm 
= constant  used  in  MTS  flame  model  (see  Eq.  2.3.2) 
= constant  used  in  MTS  flame  model  (see  Eq.  2.3.3) 


B_ 


nondimensional  depressurization  rate  coefficient 
(see  Tab.  8) 


= nondimensional  deradiation  rate  coefficient  (see 
Tab.  8) 


C 

d 

E 

f 

F 

g 

H 

i 

i. 


(0, 


(X, 


= specific  heat,  cal/g-K 

= layer  thickness,  cm 

= activation  energy,  cal/gmole 

.-0  )=  static  restoring  function  (see  Eq.  2.4.21) 
f s s 

= nondimensional  radiant  flux  intensity 

0.  -0  ) = nonautonomous  function  (see  Eq.  2.4.22) 

i,s  s ^ 

= nondimensional  surface  heat  release 


= radiant  flux  intensity,  cal/cm  -s 


o 

TTi 

n 


radiant  flux  intensity  impinging  at  the  propellant 
surface,  cal/cm^-s 


= average  molecular  weight,  g/gmole 


= exponent  in  ballistic  burning  rate  law;  also;  order 
of  the  approximating  polynomial  (see  Eq.  2.4.9) 


N. 


P 

P 


cp 

^g,s 


fl 

Q, 


r 

01 


R 

% 

s 

t 

T 


= transparency  factor  (see  Eq.  2.2.2) 

= ammonium  perchlorate  percent,  % 

= pressure,  atm 
= nondimensional  pressure 

2 

= energy  flux  intensity,  cal/cm  -s 
= nondimensional  heat  feedback  from  the  gas  phase 
= nondimensional  flame  heat  release 
= flame  heat  release,  cal/g 
= surface  heat  release,  cal/g 
= surface  reflectivity 
= universal  gas  constant,  cal/gmole-k 
= burning  rate,  cm/s 
= nondimensional  burning  rate 
= volume  scattering  coefficient,  era" 

= time,  s 
= temperature,  K 


/ 


I 
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T 

u 

ux 

u 

U 

x 

X 

w 


= nondimensional  temperature  (see  Eq.  2.2.1b) 

= nondimensional  finite  size  disturbance  of  tempe- 
rature (see  Eq.  2.4.1) 

= nondimensional  finite  size  disturbance  of  thermal 
gradient  (see  Eq,.  2.4.8) 

= gas  velocity,  cm/s 

= nondimensional  gas  velocity 

= space  variable,  cm 

= nondimensional  space  variable 

= exponent  in  KTSS  surface  pyrolysis  law  (see  Eo.2.2.3b) 


Greek  Symbols; 


* 

& 

A 

£ 

9 

\ 


f 

X 


= thermal  diffusivity,  cm  /s 
= nondimensional  layer  thickness 

= finite  difference  of  a quantity  evaluated  between 
X = 0 and  X = -|(see  Eq.  2.4.14) 

= small  quantity 

= nondimensional  temperature  (see  Eq.  2.2.1a) 

= thermal  conductivity,  cal/cm-s-K; 
also:  wavelength,  pm 

= nondimensional  thickness  of  disturbance  thermal 
layer  (see  Eq.  2.4.11) 

= density,  g/cm^ 

= nondimensional  time 

= nondimensional  time  parameter  defined  in  Eqs . 

2.32  - 2.3.4 


Subscripts  and  Superscripts: 


B 

c 

cr 

s 

f 

g 

i 

k 


= ambient  value; 

also:  absorption  layer 

- binder 

= condensed  phase  value 
= crystal 

= penetration  depth 

= final  value; 
also:  flame 

= gas 

= initial  value 
= matching  value 
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i 

= 

eneryy  loss 

n 

= 

minimum  value  for  chemical  reactions  to  occur 

P 

= 

pressure 

r 

- 

radiant 

s 

= 

surface  value 

S 

= 

static 

di 

= 

diffusion 

ki 

= 

kinetic 

re 

= 

reaction 

th 

= 

thermal 

ext 

= 

external 

ref 

= 

reference  value 

- 

= 

steady  state  or/average  value  of  a parameter 

V 

= 

vaporization 

— 00 

= 

condition  far  upstream 

+ oo 

= 

condition  far  downstream 

Abbreviation! 

5 ■ 

A 

= 

ammonia  (NH^) 

AFSC 

= 

Air  Force  System  Command 

AIAA 

= 

American  Institute  of  Aeronautics  and  Astronautics 

AIMETA 

Associazione  Italiana  di  Mecca'nica  Teorica  ed 
Applicata 

AMS 

= 

Aerospace  and  Mechanical  Sciences  Department 

AP 

= 

ammonium  perchlorate  (NH^CIO^,) 

ATI 

= 

Associazione  Termotccnica  Italiana 

BC 

= 

boundary  condition 

BC1 

= 

boundary  condition  at  x = 0 

BC2 

= 

boundary  condition  at  x = - oo 

BRL 

= 

Ballistic  Research  Laboratory 

CPIA 

a 

Chemical  Propulsion  Information  Agency 

DB 

a 

double  base 

GDF 

a 

granular  diffusion  flame  (model  of  steady  state 
burning) 

IC 

a 

initial  condition 

JSR 

= 

Journal  of  Spacecraft  and  Rockets 

KTSS 

a 

Krier-T ' ien-Sirignano-Summerf ield 

HTPB 

a 

hydroxyl- terminated  polyjautadiene 
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MTS  = Merkle-Turk-Summerf ield 

NC  = nitrocellulose 

NOTS  = Naval  Ordnance  Test  Station 

NWC  = Naval  Weapons  Center 

ODE  = ordinary  differential  equation 

ONERA  = Office  National  d'Etudec  et  de  Recherches 

Adrospatiales 

PBAA  = polybutadiene-acrylic  acid 

PA  = perchloric  acid  (HCIO^) 

PDE  = partial  differential  equation 

PMTF  = Journal  of  Applied  Mechanics  and  Technical  Physics 

PU  = polyurethane 

SRI  = Stanford  Research  Institute 

UTC  = United  Technology  Center 
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Tab. 

Tab. 

Tab. 

Tab. 

Tab. 

Tab. 

Tab. 

Tab. 


LIST  OP  TABLES 


1 Properties  of  solid  composite  propellant  AP/PBAA 
No,  941  used  as  datum  case  in  this  study. 

2 Steady  surface  temperature  dependence  on  surface 
heat  release  in  the  pressure  range  10  to  GO  atm. 

3 Surface  temperature  and  burning  rate  associated 
with  roots  A and  B,  in  the  pre'ssure  range  1 0 to 
60  atm,  for  MTS  flame.  Root  A (steady  reacting 
solution)  is  independent  on  the  order  of  the  ap- 
proximating polynomial,  while  root  B (unsteady 
reacting  solution)  does  depend.  Nonlinear  static 
restoring  function  evaluated  for  standard  condi- 
tions (T  = 300  K,  Q = -158.2  cal/g,  adiabatic 
burning) f 

4 Roots  associated  with  the  nonlinear  static  resto- 

ring function,  for  MTS  flame  with  n = 3 , at  stan- 
dard conditions  (P=30  atm,  Tg=  300  I<,  adiabatic 
burning)  showing  the  existence  of  three  static 
regimes.  For  all  cases  root  C,  corresponding  to 
the  trivial  unreacting  solution  9 = 0,  is  also 

found  and  is  always  stable. 

5 Range  of  values  of  surface  heat  release,  for  which 
self-sustained  oscillating  combustion  is  expected, 
in  function  of  pressure.  Nonlinear  static  resto- 
ring function  evaluated  for  MTS  flame  with  n=3 

at  standard  conditions  (T  = 300  K,  adiabatic  bur- 
ning) . 

6 Static  stability  limits,  in  the  pressure  range  10 
to  60  atm  for  MTS  flame  at  standard  conditions 

(r  = 300  K , Q = -158.2  cal/g),  evaluated  directly 
f?om  heat  feeSback  law. 

7 Static  stability  limits,  in  the  pressure  range  10 
to  60  atm  for  MTS  flame  at  standard  conditions 

(T  = 300  K f Q = -158.2  cal/g) t evaluated  through 
restoring  function  with  n=3. 

8a  Stabilizing  effect  of  large  radiant  flux  and  de- 
stabilizing effect  of  large  surface  heat  release 
on  stability  strength  of  steady  reacting  solution 
measured  according  to  Lyapunov.  Nonlinear  static 
restoring  function  evaluated  for  MTS  flame  with 
n=3  at  standard  conditions  (P=30  atm,  T = 300  K) . 
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Tab.  8b  Compared  stabilizing  effect  of  large  pressure 
on  stability  strength  of  steady  reacting  solu- 
tion measured  according  to  Lyapunov.  Nonlinear 
static  restoring  function  evaluated  for  MTS,  KTSS 
nonlinearized  and  KTSS  linearized  flames  with 
n = 3 at  standard  conditions  (Q  --158.2  cal/g, 

T = 300  K,  adiabatic  burning). s 

a 

i 

Tab.  9 Forcing  laws  used  for  computer  simulated  dera- 

diation and  depressurization  transients. 

Tab. 10  Computer  simulated  go/no-go  transient  tests  showing 

agreement  with  the  predicted  lower  dynamic  bounda- 
ry and  strong  dependence  of  dynamic  extinction • on 
the  initial  conditions.  All  runs  performed  accor- 
ding to  MTS  flame  for  adiabatic  opaque  strands 
with  T = 300  K and  Q = -158.2  cal/g. 


TABLE  1 


Properties  of  solid  composite  propellant  AP/PBAA 
No.  941  used  as  datum  case  in  this  study. 


ASSt'-ir.D  OR  MEASURED  PROPERTIES 

Er  ■ . crystal,  transition  heat,  Qcr 

End  . AP  vaporiz.  heat,  0V  , p 
FxuCh.  AP  decomposition  heat, 

Endoth.  binder  vapor,  heat,  Qv 
AP  cop ■ ent , p ' 

Ba  1 1 i t i c exponent , n 
KTSS  pyrolysis  law  power,  w 
Surface  activation  energy,  E 
Flame  activation  energy,  E ^ 

Condensed  phase  density,  p 
Condensed  phase  specific  heat,  C 
Condensed  phase  thermal  diffusivfty  d 
Gas  phase  specific  heat,  C 

Gas  phase  thermal  conducivlty,  X ^ 

Average  products  molecular  weight,  Ih 
Reflectivity  of  propellant  surface,  r 
Optical  absorption  coefficient,  a 
Minimum  surface  temp,  for  reactions,  T 
Matching  surface  temp,  for  pyrolysis, 

EVALUATED  PROPERTIES 

Cond. phase  thermal  conductivity,  Xc 
Refer,  frequency  factor,  A,./rGf 
Surface  gasification  heatj^Qg 
Chemical  time  constant,  Aj,j 
Diffusion  time  costant,  BM 

REFERENCE  PROPERTIES 

Pressure,  Pref 

Temperature,  T f 

Burning  Rate,  Rref  = R (pref) 

Surface  temperature,  Ts  f = t (P  ) 
Flame  temperature , Tf  s re^ 

Distance, xref  = C/Rref 
Time,  tref  = ^C/Rref 


20 
526 
800 
225 
80% 
0.46 
6.38 
1 6000 
20000 
1 .54 
0.33 
1 .4x10 
0.33 
1x10-4 
26 
O 


cal/g  of  AP 
cal/g  of  AP 
cal/g  of 
cal/g  of  binder 


cal/g mole 
cal/gmole 
g/cm3 
cal/g-K 
J cm2/s 
cal/g-K 
cal/cm-K-s 
g/gmole 


7.12x10~4  cal/cm-K-s 
2705  cm/s 

-158.2  cal/g  (+endoth.) 
0.338 
2.350 


68  atm 
300  K 

0.837  cm/s 
100O  K 
2430  K 

1 . 673x10“ 3 cm 
1 .998xlO-3  s 


Heat,  Qref  = Cc(Ts,ref  ~ Tref)  231  cal/g 

Energy  flux,  Iref  = fc  cc  Rref  (Ts,  ref^ref^  297,8  cal/cm 
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TABLE 


-126  0.863  I 0.423 


Io,  cal/cm^-s 

(X) 

Q , cal/g 

S 

(df/de  )- 

s ©. 
l.s 

0 

-158.2 

-3.53 

10 

-158.2 

•3.89 

20 

-158.2 

4.03 

30 

-158.2 

-4.27 

40 

-158.2 

-4.  50 

60 

-158.2 

-5.18 

80 

-158.2 

-6.01 

-10 

-158.2 

-3.33 

-20 

-158.2 

-3.14 

-30 

-158.2 

-2.78 

-40 

-158.2 

-2.33 

-50 

-158.2 

-1.85 

-60 

-158.2 

-1.50 

-70 

-158.2 

-1.35 

-80 

-158.2 

-0.52 

0 

-110 

-0.71 

0 

-120 

-2.10 

0 

-130 

-3.02 

0 

-140 

-3.59 

0 

-150 

-3.74 

0 

-160 

-3.51 

0 

-170 

-2.86 

0 

-180 

-1.74 

0 

-190 

0 

0 

-200 

2.41 

0 

-210 

5.85 

0 

-220 

10.60 

(x)  negative  values  represent  heat  loss  from  the  burning  surface 


TABLE 


TABLE  10 


(*)  final  pressure  in  the  case  of  depressurization  transients. 
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(MTS  flame,  n = 3) . 

- Stabilizing  effect  of  residual  radiant  flux  intensity 
on  lower  dynamic  stability  boundary  (MTS  flame,  n = 3). 

- Destabilizing  effect  of  large  surface  energy  release 
on  upper  dynamic  stability  boundary  (MTS  flame,  n = 3). 
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Fig.  15b 

Fig.  16 

Fig.  17 
Fig.  18 
Fig.  19 
Fig.  20 
Fig.  21 

Fig.  22 

Fig.  23 

Fig.  24 

Fig.  25 
Fig.  26 
Fig.  27 
Fig.  28 
Fig.  29 
Fig.  30 

Fig.  31 

Fig.  32 


Coalescence  of  A-D  roots  at  Q = -190  cal/g;  the  steady 
reacting  mode  is  stationary  for  ]Q.|  < 190  cal/g  and 
self-sustained  oscillating  for  |Q  | >190  cal/g.  (MTS 
flame,  n = 3) . See  Fig.  24. 

Negligible  influence  of  ambient  temperature  on  nonlinear 
static  restoring  function  (for  the  indicated  range  of 
values).  (MTS  flame,  n = 3) . 

Stabilizing  effect  of  pressure  on  both  reacting  modes. 
(KTSS  nonlinearized  flame,  n = 3)  . 

Effect  of  surface  heat  release  on  static  restoring  fun- 
ction according  to  KTSS  nonlinearized  flame  (n  = 3). 

KTSS  linearized  flame  (n  = 3)  applying  for  surface  tempe 
ratures  near  or  above  the  steady  reacting  value. 

Compared  static  restoring  functions  (n  = 3)  for  MTS, 

KTSS  nonlinearized  and  KTSS  linearized  flames. 

Effect  of  the  approximating  polynomial  order  on  the 
static  restoring  function  at  20  atm  of  pressure  (MTS 
flame) . 

Effect  of  the  approximating  polynomial  order  on  the 
static  restoring  function  at  30  atm  of  pressure  (MTS 
flame)  . 

Effect  of  the  approximating  polynomial  order  on  the 
static  restoring  function  at  40  atm  of  pressure  (MTS 
flame)  . 

Computer  simulated  pressurization  tests  confirming  the 
existence  of  three  possible  static  regimes  (stationary 
for  Qu  = -180  cal/g,  self-sustained  oscillating  for 
Q = - 200  cal/g,  extinguished  for  Qs  = - 220  cal/g) . 

Cf.  Tab.  4. 

Static  instability  of  burning  propellants  for  large 
surface  heat  release. 

Approximate  construction  of  the  static  stability 
boundary  and  its  meaning  (MTS  flame) . 

Rigorous  construction  of  the  static  stability  boundary 
for  Qs  = - 120  cal/g  (MTS  flame). 

Rigorous  construction  of  the  static  stability  boundary 
for  O = - 158.  2 cal/g  (MTS  flame). 

Rigorous  construction  of  the  static  stability  boundary 
for  0S  = - 200  cal/g  (MTS  flame) . 

Evaluating  the  static  stability  limit  through  the 
static  restoring  function  at  20  atm  of  pressure  (MTS 
flame , n = 3) . 

Evaluating  the  static  stability  limit  through  the 
static  restoring  function  at  30  atm  of  pressure  (MTS 
flame,  n = 3) . 

Evaluating  the  static  stability  limit  through  the 
static  restoring  function  at  40  atm  of  pressure  (MTS 
flame,  n = 3) . 
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Fig.  33a  - Representative  time  histories  of  burning  rate  during 
deradiation  showing  possible  occurence  of  dynamic 
extinction. 

Fig.  33b  - Corresponding  trajectories  in  burning  rate  vs  heat 
feedback  plane. 


Fig.  33c 
Fig.  33d 


- Nomenclature  used  for  deradiation  transients. 


Nature  of  nonautonomous  function  considered  in  this 
study . 

Fig.  34  - The  range  of  possible  unstable  equilibrium  points  in' 

dynamic  burning  regime  is  limited  by  the  no-return 
point  M (the  unstable  root  associated  to  the  static 
• restoring  function  for  *-*•<*>). 

- Static  and  dynamic  stability  boundaries  on  a burning 
rate  vs  heat  feedback  plot  (MTS  flame) . 

- Static  and  dynamic  (MTS  flame,  n = 3)  stability 
boundaries  on  a burning  rate  is  pressure  plot. 

- Go/no-go  computed  deradiation  tests  showing  occurrence 
of  dynamic  extinction  (Br  = 5) . 

- Go/no-go  computed  deradiation  tests  showing  occurrence 
of  dynamic  extinction  (Br  = 10) . 

- Dynamic  extinction  following  linear  deradiation. 

- Dynamic  extinction  following  bilinear  deradiation. 

- Dynamic  extinction  following  parabolic  deradiation. 

- Dynamic  extinction  following  exponential  deradiation 
at  10  atm  of  pressure. 

- Dynamic  extinction  following  exponential  deradiation 
at  20  atm  of  pressure. 

- Dynamic  extinction  following  exponential  deradiation 
at  30  atm  of  pressure. 


Fig.  75 

Fig.  36 

Fig.  37 

Fig.  38 

Fig.  39 
Fig.  40 
Fig.  41 
Fig.  42 

Fig.  43 

Fig.  44 

Fig.  45 

Fig.  46 

Fig.  47 

Fig.  48a 


- Dynamic  extinction  following  exponeptrial  deradiation 
at  40  atm  of  pressure. 

- Dynamic  extinction  followincp'^xponential  depressurization 
from  a par&metically  changed  initial  pressure. 

- Dynamic  extinction  f^ilowing  exponential  depressurization 
with  a parametrically  changed  depressurization  rate. 

- Sketch  of  tho^shock  tube  apparatus. 


Fig.  48b  - Sketch^crf  the  shock  tube  test  section. 

Fig.  49  - SJfesftic  pressure  deflagration  limits  of  the  composite  pro- 

^pellant  as  measured  in  the  shock  tube  test  section. 


Fig^-SO 


- Composite  propellant  tested  in  the  shock  tube.  Ignition 
at  1 atm  of  nitrogen.  UDper  trace  : pressure  (10  atm/div) , 
lower  trace  : emitted  radiation  (relative  intensity) . 

Time  scale  : 5 ms/div. 
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Fig.  51  - See  above.  Time  scale:  1 ms/div. 

Fig.  52  - Noncatalyzed  DB  tested  in  the  niston  tube.  Ignition  at 
0.50  atm  of  air.  Upper  trace  : pressure  (20  atm/div), 
lower  trace  : emitted  radiation  (relative  intensity) . 
Time  scale  : 5 ms/div. 

Fig.  53  - See  above.  Time  scale  : 50  ms/div. 

Fig.  54  - Noncatalyzed  DB  tested  in  the  piston  tube.  Ignition  at 
0.75  atm  of  air.  Upper  trace  : pressure  (20  atm/div), 

• lower  trace:  emitted  radiation  (relative  intensity) . 

Time  scale  : 5 ms/div. 


Fig . 


55 


See  above.  Time  scale 


50  ms/div. 
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Fig.  33a  Representative  time  histories  of  burning  rate  during  deradi- 
ation showing  possible  occurrence  of  dynamic  extinction. 

Fig.  33b  Corresponding  trajectories  in  burning  rate  vs  heat  feedback  plane 

Fig.  33c  Nomenclature  used  for  deradiation  transients. 
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CHAPTER  3 

APPLICATION  OF  THE  LASER  DOPPLER  VELOCI  METRY  TO 
THE  FLAME  STRUCTURE  OF  SOLID  PROPELLANTS 

,Sec.  3.1  - INTRODUCTION 

Sec.  3.2  - LASER  DOPPLER  VELOCIMETRY 

Sec.  3.3  - MIE  SCATTERING 

Sec.  3.4  - EXPERIMENTAL  SET-UP 

Sec.  3.5  - EXPERIMENTAL  RESULTS 

Sec.  3.6  - CONCLUSIONS  AND  FUTURE  WORK. 

References 
Table  1 

List  of  figures 
Figures  1-7 
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Sec.  3.1  - INTRODUCTION 

The  LDV  technique  is  intaided  to  measure  the  velocity 
profiles  in  the  gaseous  region  near  the  burning  surface  of  a 
solid  propellant.  The  main  reason  of  this  choice  is  the  fact 
that  the  usual  intrusive  experimental  techniques  are  of  lit- 
tle help  in  a high  temperature  reactive  medium  with  large 
thermal  gradients.  By  other  standard  techniques  the  flame  tem- 
perature, the  burning  rate  and  the  condensed  phase  thermal 
profile  can  also  be  determined.  This  would  allow  to  cross- 
check the  LDV  results  in  the  gas  phase  zone. 

A steady  state  strand  burner,  with  two  symmetrical  and 
opposite  optical  windows,  was  specifically  designed  for  LDV 
experiments.  The  operating  pressure  range  is  1 to  1 5 atm.  A 
computerized  data  acquisition  and  processing  system  was  also 
realized.  The  facility  has  been  applied  to  a steadily  burning 
solid  propellant,  with  the  purpose  of  exploring  the  possibility 
of  future  applications  to  unsteady  situations. 

Several  tests  of  preliminary  nature  were  performed  using 
double-base  (DB)  (both  catalized  and  noncatalized)  and  ammo- 
nium perclorate  (AP)  based  composite  solid  propellants.  The 
differential  mode  of  operation  was  used  for  the  LDV  instru- 
ment, with  observation  in  the  forward  direction.  The  particles 
present  in  the  plume  of  the  burning  solid  propellant  were  used 
as  scattering  centers.  Reasonably  good  results  of  the  LDV  mea- 
surements were  obtained  in  the  steady  burning  of  composite  pro- 
pellants; in  principle,  useful  information  could  be  obtained 
in  unsteady  situations  as  well.  Some  problems  remain  to  be 
solved  and  more  data  are  needed  in  order  to  cross-check  the 
LDV  results  and  to  clarify  the  uncertainties  so  far  encounte- 
red. 


129 


Sec.  3.2  - LASER  DOPPLER  VELOCIMETRY 

* The  basic  principle  involved  in  Laser  Doppler  Velocimetry 

(LDV)  is  the  Doppler  frequency  shift  of  a monochromatic  and 
coherent  light  beam  scattered  by  micrometric  particles  suspen- 
ded in  a moving  medium. 

The  Doppler  frequency  shift  in  normally  resolved  by  hete- 
rodyning the  scattered  light  with  a reference  beam  on  the 
surface  of  a photodetector  (reference  beam  LDV) . Alternative- 
ly, scattered  waves  coming  from  two  incident  beams  can  be 
heterodyned  (differential  LDV) . The  resulting  light  intensity 
and,  hence,  the  photocurrent  will  be  modulated  at  the  Doppler 
frequency  (Refs.  1-3). 

(3.2.1)  f = 21 J S i M ( p/2)  / \ 

A P 


where  U is  the  particle  velocity  component  in  the  direction 
of  the  bisector  of  the  illuminating  cross-beams;  k = 2 tf/X 
is  the  wave  vector;  X is  the  wave-length  of  the  laser  beam  and 
J3  is  the  angle  between  the  two  cross-beams  (see  Fig.  1). 

LDV  technique  allows  a direct  measurement  of  velocity 
components  (through  the  measurement  of  the  Doppler  frequency 
by  an  electronic  signal  processor)  with  a very  sharp  space 
resolution  (tfO.1  mm3)  an(j  without  disturbing  the  flow.  This 
is  an  important  feature  in  reactive  media.  However,  asuffi- 
ciently  large  number  of  scattering  particles  (^10?  M ) is 
required  to  obtain  continous  velocity  information. 

The  performance  of  an  LDV  system  strongly  depends  on 
the  size  distribution  and  concentration  of  the  particles  scat 
tering  light  into  the  photodetector.  Particle  concentration 
has  to  be  high  to  obtain  good  time  resolution  of  the  velocity 
variation,  while  particle  sizes  have  to  be  very  small  to 
follow  spatial  velocity  gradients.  The  capability  of  particles 
to  follow  the  gas  flow  can  be  roughly  estimated  from  the 
Bassett's  general  equation  (Refs.  4 and  5).  In  unsteady  flows 
the  LDV  chain  of  instruments  must  be  treated  as  a transducer 
for  which  the  transfer  function  is  determined  by  the  dynamics 
of  tracer  particles. 

Parallel  measurements  of  particle  size  and  density  are 
of  the  greatest  importance  in  experiments  where  flow  seeding 
with  particles  of  known  size  is  not  feasible  or  it  cannot  be 
assumed  "a  priori"  that  ambient  particles  are  adeguately 
following  the  fluid  flow,  LDV  systems  can  provide  informations 
on  particle  size  from  the  correlation  of  particle  diameter 
with  the  shape  of  the  LDV  signal. 


Sec.  3.3  - MIE  SCATTERING 

It  was  shown  (Refs.  6-9)  that  the  single  particle  Doppler 
signal  can  be  accurately  predicted  on  the  basis  of  the  Mie 
scattering  theory.  The  general  expression  of  the  Doppler 
photocurrent,  produced  by  a single  particle  crossing  the  geo- 
metric center  of  the  probe  volume,  is  given  by 


(3.3.1) 


k 


where  ^ is  the  detector  sensitivity,  i>*D  = 2 7T  f^,  and  the  two 
j r -i lent  beams  are  assumed  to  have  equal  intensities 
If'  " I02  = 1<y  T*le  tcrm  p represents  the  "pedestal"  amplitude 
oi'the  current  (Fig.  2),  while  D represents  the  Doppler  ampli- 
tude and  ^ determines  the  phase  of  the  scattered  intensity.  P, 

D and  are  integral  quantities,  integrated  over  the  collecting 
solid  angle, -fl  , of  the  receiving  optics. 

They  depend  an  the  direction  of  detection  (<^,  </  ),  i.e.  the  axis 
of  the  collecting  aperture,  end  are  defined  in  terms  of  the 
complex  amplitude  functions  given  by  the  Mie  theory  (Refs.  10  - 
12). 

The  quantities  P,  D and  -'fr  depend  only  on  the  scattering 
properties  and  size  of  the  particle  (for  a fixed  LDV  geometry) , 
whereas  thenhotocurrent  is  a time-varying  function  of  the 
particle  position  in  the  probe  volume.  Indeed,  a moving  fringe 
pattern  is  seen  by  the  detector  due  to  two  beam  interference 
and  particle  motion. 


The  fringe  contrast  or  visibility  is  defined  by 


(3.3.2) 


V = D/P 


The  exact  dependence  of  V on  particle  diameter,  d,  can  be  nu- 
merically obtained  through  P and  D evaluated  by  a computer  code 
(Ref.  9)  based  on  Mie  scattering  formula.  Visibility  is  an 
important  parameter  because  its  computed  values  can  be  easily 
compared  with  experimental  ones. 

Experimental  results  reported  in  the  literature  showed 
the  possibility  of  parallel  measurements  of  individual  particle 
velocity  and  size  (Refs.  7,  8,  9 and  13).  Applications  include 
particle  size  measurements  from  calibrated  monodisperse  aerosols, 
polydisperse  fuel  sprays  in  combustion  and  dusts  in  the  atmo- 
sphere. Comparisons  with  other  particle  sizing  methods  indicate 
that,  within  the  bounds  of  experimental  error,  the  use  of  the 
visibility  function  yields  particle  size  distribution  measure- 
ments which  agree:  with  the  results  obtained  by  other  methods. 

Unfortunately,  the  visibility  signal  cannot  be  related  to 
particle  size  over  a wide  range  by  a monotonic  function  (see 
Fig. 3 ).  Therefore,  some  kind  of  correlation  of  particle  visi- 
bility with  signal  magnitude  or  pedestal,  P,  is  also  necessary. 
Moreover,  the  visibility  is  a sensitive  function  of  the  index 
of  refraction  of  particles.  If  this  parameter  is  unknown,  some 
other  technique  has  to  be  used.  For  example  the  measurements  of 
the  ratio  of  the  light  scattered  from  individual  particles  at* 


two  or  more  finite  angles.  Combined  results  of  the  two  techni- 
ques have  been  proved  useful  in  making  "in  situ"  measurements 
of  particle  diameters. 
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Sec.  3.4  - EXPERIMENTAL  SET-UP 

The  LDV  measurements  were  performed  in  a steady  state 
strand  burner,  with  two  symmetrical  and  opposite  optical 
windows,  specifically  designed  for  this  purpose.  In  the 
final  version  of  the  experimental  set-up  (see  Fig.  4),  two 
thick  lenses  were  mounted  directly  on  the  burner.  Advanta- 
ges of  this  configuration  are  simplicity  and  the  possibili- 
ty of  utilising  all  of  the  small  window  aperture  that  is 
available . 

The  differential  LDV  mode  was  chosen  because  of  its 
higher  signal  to  noise  characteristics  at  moderate  particle 
concentrations.  Moreover,  this  type  of  LDV  system  is  the 
best  suited  for  individual  realization  velocimetry  and  paral- 
lel particle  size  analysis.  It  is  not  difficult  to  realize 
the  situation  in  which  the  Doppler  signal,  that  is  available 
for  processing,  is  produced  by  no  more  than  one  particle  in 
the  probe  volume  at  a time.  In  fact,  it  is  possible  to  adjust 
the  probe  volume  dimensions,  both  in  diameter  and  in  length, 
by  selecting  the  proper  cross-beam  angle  p and  the  magnifi- 
cation of  the  light  collecting  system.  Typical  values  relati- 
ve to  our  experiments  are  given  in  Table  1 . The  main  advanta- 
ge of  the  differential  LDV  is  that  it  is  quite  simple  to  align 
and  is  not  sensitive  to  small  vibrations.  Also,  the  Doppler 
frequency  is  independent  of  the  detection  angle  and  a large 
collection  aperture  can  be  used  without  a spread  of  Doppler 
frequency . 

The  LDV  system  comprises  a 5 mW  He-Ne  laser  and  a beam 
splitter  with  variable  beam  separation,  which  allows  conti- 
nuous variation  of  the  cross-beam  angle.  The  two  inci- 

dent parallel  laser  beams  were  directed  toward  the  lens  L^ , 
in  such  a way  that  one  laser  beam  was  centered  on  the  lens 
axis.  This  beam  is  not  diverted  from  its  horizontal  direction, 
but  focused  on  the  focal  point  (on  the  axis  of  the  strand  bur- 
ner) . The  second  beam  is  impinging  on  the  lens  L^  at  a distan- 
ce h from  the  axis,  and  is  deflected  and  focused  in  the  same 
focal  point.  The  crossing  region  defines  the  probe  volume  of 
the  LDV  system  and  the  velocitv  component  in  the  vertical  di- 
rection can  be  measured,  With  the  geometry  of 

Fig.  4 the  measured  velocity  component  is  not  exactly  the  axial 
one;  however  the  difference,  of  the  order  of  sin  {p/2),  is 
not  significant  for  small  p. 

LDV  experiments  were  performed  by  using  as  scattering  cen- 
ters the  particles  present  in  the  plume  of  the  burning  solid 
propellant  sample. No  external  seeding  was  used  in  order  to 
avoid  any  disturbance  of  the  combustion  processes.  Prelimina- 
ry experiments  showed  that,  generally,  sufficiently  high  parti- 
cle rates  are  found  in  the  region  near  the  burning  surface. 

The  experimental  conditions  are  now  schematically  descri- 
bed. Initially,  the  rod  of  solid  propellant,  centered  on  the 
axis  of  the  burner,  shuts  out  the  two  incident  laser  beams. 
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When  the  rod  burns,  its  surface  goes  down  and  at  the  instant 
tc  the  horizontal  beam  is  allowed  to  pass.  After  a delay  At, 
the  second  beam  is  also  passing  and  the  scattered  light  with 
Doppler  informations  is  received  by  the  photomultiplier . A 
triggering  signal,  coming  out  from  the  photodiode  (P.D.)  in- 
tercepting the  horizontal  beam,  precisely  defines  the  instant 
t0 . From  this  the  axial  position  of  the  LDV  probe  volume,  with 
respect  to  the  burning  surface,  can  be  inferred  by  comparison 
of  the  time  delay  and  the  propellant  burning  rate.  The  strand 
burner  is  also  equipped  for  simultaneous  measurements  of  the 
. propellant  burning  rate  and  the  instantaneous  pressure  (Ref. 

14)  . 

The  receiving  optics  of  the  LDV  system  comprises:  the  lens 
, fixed  on  the  burner  window,  that  collects  the  radiation  scat- 
tered in  the  forward  direction;  the  lens  L,  that  is  movable  and 
focuses  the  scattered  light  on  a 0.3  mm  pinhole  in  front  of  a 
photomultiplier.  This  is  equipped  with  an  interferential  opti- 
cal filter,  centered  at  632.8  +1.0  nm  (the  laser  wave-length), 
in  order  to  reject  flame  emission. 

The  photomultiplier  signal  is  processed  by  an  electronic 
counter  processor  (DISA  mod.  55L90)  that  allows  working  with 
variable  particle  concentrations  and  does  not  have  dropout  pro- 
blems. Moreover,  it  has  a large  dynamic  range  and  no  slewrate 
limitations;  it  accepts  individual  signals  and  determine  the 
correct  Doppler  frequency,  hence,  the  particle  velocity.  Pro- 
per use  of  this  instrument  would  require  rejection  of  multiple 
particle  signals^,  because  of  random  phase  fluctuations  which 
will  lead  to  incorrect  velocity  measurements  (Ref.  15).  This 
requirement  can  be  generally  satisfied  by  proper  reduction  of 
the  probe  volume  dimensions. 


134 


Sec.  3.5  - EXPERIMENTAL  RESULTS 

A typical  individual  Doppler  signal,  after  amplification 
and  band-pass  filtering,  is  shown  in  Fig.  5.  Many  of  these 
signals  were  recovered  in  order  to  optimize  the  optical  sy- 
stem and  to  define  the  effective  probe  volume  dimension, 
by  comparing  the  observed  and  the  computed  fringe  number. 

With  the  experimental  conditions  summarized  in  Tab.  1 
and  the  solid  propellants  used,  we  observed  that  generally 
no  more  than  one  particle  at  once  was  present  in  the  probe 
volume.  If  the  pass-band  filtered  Doppler  signal,  relative 
to  each  particle,  satisfies  a number  of  validation  conditions 
(in  terms  of  amplitude  level  and  signal  to  noise  ratio) , the 
counter  processor  measures  the  Doppler  frequency  by  means  of 
a zero-crossing  method.  The  output  is  velocity  information 
that  is  retained  until  another  particle  gives  a validated  si- 
gnal . 

During  preliminary  experiments,  the  velocity  time  history 
was  monitored  by  a storage  oscilloscope  (Tektronix,  mod. 466), 
after  digital  to  analog  conversion  of  the  output  of  the  Dop- 
pler counter  processor  . Typical  results  obtained 

with  a noncatalyzed  double-base  (a,b)  and  an  AP  based  composi- 
te propellants  (c,  d)  are  reported  in  Fig.  6.  The  upper  trace 
is  the  continuous  recording  of  the  photodiode  output;  the  lo- 
wer trace  is  the  velocity  recording  in  which  each  point  corre- 
sponds to  a single  particle  measurement.  Oscilloscope  records 
are  triggered  by  the  photodiode  signal  and  have  the  same  time 
base.  It  can  be  noted  that  large  variations  with  sharp  gradients 
characterise  the  velocity  of  the  DB  propellant  ir  conjunction 
with  irregular  fluctuations  in  the  particle  rate  (i.e.  the 
number  of  points/div.  in  velocity  traces)'.  At  exactly  the  same 
experimental  conditions,  higher  particle  rates  and  more  regular 
velocity  traces  were  found  with  the  AP  propellant. 

The  time  resolution  of  the  LDV  measurements  depends  on  the 
particle  rate,  n,  hence  on  the  particle  concentration  N,  the 
mean  particle  axial  velocity  U and  the  probe  volume  cross-sec- 
tion S: 


(3.5.1) 


n = cl  N US 


where  <*  is  an  efficiency  factor  taking  into  account  the  valida- 
tion rate  of  particle  signals  from  the  counter  processor.  In 
our ^experiments  with  double-base  propellants  the  time  history 
of  n was  measured  by  a multiscalear.  Many  tests  were  made  during 
preliminary  velocity  measurements  and  the  mean  values  of  n for 
different  optical  and  electronic  configurations  of  the  LDV  were 
determined  in  order  to  match  the  optical  system  to  the  mean 
particle  concentration.  Large  fluctuations  in  n were  found  from 
0 to  a maximum  value  of  about  1 Ktfzjthat  will  correspond  to  a 
particle  concentration  of  about  2x1 08  m“3. 
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With  these  experimental  conditions  the  LDV  cannot  provide 
a continuous  record  of  the  gas  velocity.  Moreover,  since  the 
particles  cross  randomly  the  probe  volume,  they  provide  a ran- 
domly sampled  time  series  of  the  fluctuating  velocity.  Time  re- 
solution can  be  limited  by  this  characteristic.  This  is  an  im- 
portant feature  for  unsteady  situations  and  rapidly  fluctuating 
flows.  As  far  as  a steady  state  situation  is  concerned,  the 
particle  rate  is  more  than  adequate  to  velocity  measurements 
by  LDV.  Present  results  are  limited  to  this  case. 

Data  reduction  was  performed  by  a minicomputer  (WANG,  mod. 
2200) . The  computer  interface,  which  reads  the  velocity  data 
supplied  by  the  counter  processor,  was  activated  by  a trigge- 
ring pulse  derived  by  the  photodiode  signal.  Once  the  data  rela- 
tive to  one  test  (about  2*000  velocity  measurements)  have  been 
acquired  into  memory,  they  are  transferred  to  a tape  for  perma- 
nent storage.  To  date,  with  the  minicomputer  at  our  disposal, 
the  maximum  data  acquisition  rate  is  of  about  2 KHz.  Analysis 
programs  have  been  developed  for  the  LDV  measurements  and  are 
described  in  Ref.  16.  An  example  of  the  results  obtained  v/ith 
this  computerized  system  is  shown  in  Fig.  7.  Points  represent 
the  mean  velocity  evaluated  in  a fixed  period  of  time  (=12.5  ms) 
over  about  10-20  particles.  The  total  time  base  covers  more 
than  1 s,  corresponding  to  the  burning  of  about  5 mm  of  the  so- 
lid propellant  sample.  This  velocity  curve  refers  to  the  same 
experimental  test  of  Fig.  6 d. 


Sec.  3.6  - CONCLUSIONS  AND  FUTURE  WORK 


Data  collected  up  to  now  suggest  the  possibility  of  re- 
sonable  measurements  of  particle  velocity  in  the  gaseous 
region  above  the  burning  surface.  So  far,  only  tests  in  a 
steady  state  situation  at  pressures  up  to  6 atm  were  perfor- 
med, but  futher  progress  is  being  made.  The  possibility  of 
instantaneous  measurements  of  gas  velocity  versus  time  is 
strictly  dependent  on  scattering  particle  concentration  and 
particle  size  distribution.  Generally,  particle  concentra- 
tion is  not  uniform  in  time  and  depends  on  the  solid  propel- 
lent type;  particle  size  distribution  is,  up  to  now,  an  unk- 
nown parameter  and  will  be  of  fundamental  importance  to  deter- 
mine it  in  order  to  be  sure  that  particles  follow  the  gas 
flow. 


The  minimum  time  resolution  that  can  be  obtained  depends 
on  the  mean  particle  velocity  (e;  15  m/s),  the  probe  volume 
dimension  in  the  axial  direction  (=0.15  mm)  and  the  particle 
rate.  If  the  typical  transit  time  for  a particle  to  cross  the 
probe  volume  is  of  about  10/as,  the  maximum  particle  rate 
could  be  100  KHz,  still  satisfying  the  condition  of  single 
particle  scattering.  This  theoretical  result  is  well  above 
the  observed  one  and,  in  our  situation,  is  quite  impossible 
to  define  a fixed  time  resolution  because  of  the  large  fluc- 
tuations of  particle  rate. 

Several  tens  of  runs  with  the  LDV  technique  show  more 
complicated  profiles  than  expected.  Although  results  are  not 
always  reproducible,  we  may  affirm  that  the  velocity  beha- 
viour presents  wide  irregular  fluctuations  in  the  case  of 
the  DB  propellants  (Fig.  6 a and  b)  and  a surprising  minimum 
after  about  200  ms  in  some  tests  with  the  AP  propellant 
(Fig.  7) . In  all  runs  a large  velocity  dispersion  (up  to  20%) 
was  observed  around  the  mean  value;  this  is  difficult  to  ex- 
plain on  the  basis  of  fluidynamics  considerations.  Compari- 
son with  shadograph  pictures  at  similar  test  conditions  sug- 
gest a quasi-monodimensional  laminar  gas  flow. 

The  discrepancy  with  the  presumed  steady  state  laminar 
flame  conditions  could  be  due  to  a polydisperse  particle  si- 
ze distribution  and  the  consequent  difficulty  for  larger 
particles  or  agglomerates  to  follow  accurately  the  gas  flow. 
Indeed,  the  measured  mean  velocities,  for  both  the  DB  and  AP 
propellant,  seem  to  be  below  (about  30%)  the  values  obtaineei. 
by  the  measured  burning  rate  of  the  condensed  phase  and  the 
estimated  flame  temperature,  assuming  a ouasi-steady  state 
mass  balance.  In  addition,  overlapping  effects  are  due  to  so- 
mewhat time-varying  pressure  during  the  burning  of  the  solid 
propellant. 

The  precision  of  the  LDV  instrument  itself  could  be  que- 
stioned, in  this  difficult  application,  but  comparative  tests 
with  cold  flows  in  the  same  burner  and  the  analysis  of  a lar- 
ge number  of  Doppler  signals  during  burning  conditions  con- 
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firmed  the  reliability  of  LDV  measurements  with  experimental 
errors  in  the  bounds  of  + 5%. 

In  conclusion,  a polydisperse  particle  size  distribution 
could  be  the  major  source  of  the  wide  spread  in  the  measured 
instantaneous  velocities.  Indeed,  the  particle  size  evaluation 
by  means  of  electron  microscope  analysis  seems  to  confirm  this 
hypothesis,  showing  a large  number  of  particle  diameter  well 
above  10 /Am.  However  the  particle  capture  could  cause  agglome- 
ration and  clustering.  The  same  phenomena  may  also  occur  whi- 
le particles  are  traveling  and  reacting  in  the  gaseous  stream 
and  this  could  explain  the  behaviour  of  Fig.  7. 

To  resolve  this  problem,  a research  program  is  starting 
with  the  purpose  of  a systematic  analysis  of  particle  size. 

The  first  attempt  will  be  performed  by  measuring  relative  scat 
tered  light  intensities  at  different  scattered  angles  and  com- 
paring with  theoretical  values  computed  on  the  basis  of  the 
Mie  theory.  Qualitative  informations  could  be  obtained,  but 
is  should  be  noted  that  this  technique  is  only  applicable  to 
particle  diameters  in  the  range  of  0.5 /Am  to  5 /4m  (Ref.i?). 
Particle  sizing  by  visibility  measurements  will  be  used  to 
measure  diameters  in  the  range  5/<m  to  50  pm.  The  second  tec- 
nique  has  been  proven  practical  in  making  in  situ  measurements 
of  particle  diameters  simultaneously  with  measurements  of  ve- 
locity and  particle  concentration  (Ref. 13)  . The  ideal  situa- 

tion will  be  characterised  by  the  possibility  of  parallel  mea- 
surement of  size  and  velocity  for  each  individual  particle. 
This  is  possible,  in  principle,  but  requires  a more  sophisti- 
cated and  expensive  electronic  instrumentation. 
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LIST  OF  FIGURES 

* 

Figure  1 : Geometry  of  a differential  LDV  system:  ^8  is  the 

cross-angle,  IT.,  and  K*.  2 the  directions  of  the 
two  incident  beams. 

Figure  2:  Doppler  signal  generated  by  a single  moving 

particle . 

Figure  3:  Typical  signal  visibility  curve  versus  particle 

diameter . 

Figure  4:  Experimental  set-up. 

Figure  5:  A typical  individual  Doppler  signal  after  ampli- 

fication and  band-pass  filtering. 

Figure  6:  Particle  velocity  behaviour  versus  time,  a)  and 

b) : DB  propellants;  c)  and  d) : AP  propellants. 

Burning  conditions:  p = 6 atm 
Scales:  t = 10  ms/div;  U=  5 m/s/div. 

Particle  velocity  versus  time  obtained  by  the 
computerized  data  acquisition  system.  Each 
point  represents  the  mean  velocity  evaluated 
in  a fixed  period  of  time  ( — 12.5  ms)  over 

about  10-20  particles*  AP  propellants. 


Figure  7: 


L 


PARTICLE 

TRAJECTORY 


X 


FIG.  5 


i 


CONCLUSIONS  AND  FUTURE  WORK 

Sec.  4.1  - CONCLUSIONS 
Sec.  4.2  - FUTURE  WORK. 


r 

148 


[ 

I i 


Sec.  4.1  - CONCLUSIONS  \ 

A nonlinear  stability  analysis  of  solid  propellant  bur- 
ning was  carried  out,  within  the  framework  of  a thermal  theo- 
ry and  for  quasi-steady  gas  chase,  allowing  for  finite  size 
disturbances.  This  required  an  integral  method  in  reducing 
the  partial  differential  equation  for  the  condensed  phase 
heat  feedback  to  an  approximate  ordinary  differential  equa- 
tion. It  is  shown  that  a nonlinear  algebraic  function,  cal- 
led restoring  function,  can  be  defined  that  contains  all  ba- 
ric properties  of  equilibrium  and  stability  of  burning  solid 
oropellants . This  function  does  not  depend  on  time,  but  only 
on  the  nature  of  the  solid  propellant  (including  its  flame) 
and  the  operating  conditions  (pressure,  ambient  temperature, 
and  energy  exchange  with  surroundings) . Analysis  of  the  non- 
linear algebraic  restoring  function  reveals  that  two  well 
defined  burning  regimes  exist,  each  limited  by  stability 
boundaries:  the  static  and  the  dynamic  regimes.  Of  these  two, 
the  region  of  dynamic  burning  is  wider,  in  that  under  dynamic 
conditions  a propellant  may  pass  through  a region  which  is 
statically  unstable  but  dynamically  stable. 

The  static  regime  can  be  observed  experimentally  and  the- 
refore can  be  studied  also  in  the  framework  of  Zeldovich  ap- 
proach. Static  stability  boundaries  are  defined  as  the  lines 
of  separation  between  regions  where  stable  steady  state  solu-  y 

tions  are  allowed  and  regions  where  only  unstable  or  no  stea- 
dy solutions  are  found.  For  usual  prooellants,  only  a lower 
(burning  rate  of  the  order  of  0.1  cm/s  or  less)  static  stabi-  i 
lity  boundary  is  found.  For  propellants  with  large  surface 
heat  release,  lower  and  UDoer  (burning  rate  of  the  order  of 
1 cm/s  or  more)  static  stability  boundaries  are  found.  Between 
these  two  static  stability  boundaries,  the  stable  steady  solu- 
tion is  either  stationary  or  self-sustained  oscillating.  Seve- 
ral methods  are  suggested  as  to  the  determination  of  the  sta- 
tic stability  boundaries. 

The  dynamic  regime  cannot  be  observed  experimentally  in 
a stationary  mode  and  can  only  be  studied  in  the  framework  of 
a flame  model.  The  lower  dynamic  stability  boundary  is  defined 
as  that  ultimate  burning  condition  beyond  which  extinction  ne- 
cessarily follows  during  a burning  transient.  The  upper  dyna- 
mic stability  boundary  implies  vigorous  accelerations  of  the 
combustion  wave,  possibly  followed  by  dynamic  extinction.  For 
usual  propellants  only  the  former  is  found;  for  propellants 
with  large  surface  heat  release  the  latter  is  also  found. 

It  is  shown  that  the  lower  dynamic  stability  boundary 
holds  true  both  for  deradiation  and  depressurization,  for  o- 
paque  as  well  as  transparent  condensed  phase,  for  fast  dece- 
leration of  the  combustion  wave  (e.g. , by  depressurization) 
as  well  as  for  fast  acceleration  (e.g.,  by  pressurization)  if 
an  excessively  large  burning  rate  overshoot  is  attained.  The 
lower  dynamic  stability  boundary  was  determined  as  an  asympto- 
tic (in  time)  boundary  for  monotonically  decreasing  external 
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controlling  parameters.  If  no  change  in  time  of  the  external 
controlling  parameters  occurs,  the  propellant  is  only  subjec- 
ted to  random  intrinsic  disturbances  and  the  static  stabili- 
ty analysis  apply.  Moreover,  if  the  effect  of  the  time  chan- 
ge of  the  external  controlling  parameters  (nonautonomous  func- 
tion) is  negligible  compared  to  the  restoring  function,  the 
lower  dynamic  stability  boundary  holds  true  for  any  time  (e- 
ven  finite)  and  for  any  external  law  (even  non  monotonic  or 
increasing) . 

MTS,  KTSS  nonlinearized  and  KTSS  linearized  unsteady  fla- 
me models  were  implemented  in  this  study.  While  KTSS  lineari- 
zed is  of  no  value  for  burning  rate  less  than  about  90S  of  the 
steady  state  value,  very  reasonable  and  similar  results  are 
displayed  by  MTS  and  KTSS  nonlinearized.  However,  MTS  flame 
model  is  considered  superior  since  it  accounts  also  for  chemi- 
cal kinetics  (but  it  requires  two  constants  to  be  evaluated) . 

The  order  of  the  polynomial  chosen  to  approximate  the  di- 
sturbance thermal  profile  in  the  condensed  phase  does  not  af- 
fect qualitatively  the  shape  of  the  restoring  function.  A cu- 
bic law  was  found  to  give  accurate  predictions  for  pressures 
up  to  about  30  atm;  a quadratic  law  seems  more  appropriate 
for  larger  pressures. 

The  validity  of  this  nonlinear  stability  theory  was  veri- 
fied by  computer  simulated  transients.  In  general,  excellent 
agreement  was  found  between  the  analytical  predictions  and 
the  numerical  results  obtained  by  integration  of  the  governing 
partial  differential  equation.  Numerical  values^yre  given  only 
for  a particular,  ammonium  percholorate  based,  composite  pro- 
pellant. It  is.fclt,  however,  that  all  analyses  were  conduc- 
ted from  a broad  point  cf  view  and,  therefore,  the  conceptual 
findings  may  be  extended  to  other  types  of  solid  propellants, 
provided  a proper  flame  model  is  employed. 

In  conclusion,  nonlinear  propellant  burning  instability 
can  be  immediately  defined  from  the  knowledge  of  the  associa- 
ted restoring  function.  This  a property  strictly  dependent  on 
the  nature  of  the  propellant  and  the  operating  conditions. 

Experimental  work  is  in  progress. 'Three  different  solid 
propellants  are  being  characterized;  an  ammonium  percholorate 
based  composite,  a catalyzed  double  - base,  and  a noncatalyzed 
double  - base.  Steady  state  burning  rates  were  measured  in  a 
strand  burner;  thermal  profiles  in  the  condensed  phase  and  ra- 
diative emissions  are  being  detected;  high  speed  movies  are 
being  taken.  Depressurization  and  pressurizations  tests  are 
being  conducted  respectively  in  a specifically  designed  strand 
burner  and  in  a properly  modified  shock  tube.  Laser  doppler  ane- 
mometry  is  being  applied  to  burning  propellants  in  order  to  be 
able  to  measure  the  gas  velocity  in  the  plume  of  the  sample. 


Sec.  4.2  - FUTURE  WORK 


Other  unsteady  flame  models  will  be  investigated,  in  par- 
ticular the  Kooker-Zinn  model  proposed  in  1974.  Results  simi- 
lar to  those  obtained  with  the  KTSS  linearized  flame  model 
are  expected.  Further  improvements  of  the  analytical  techni- 
que used  to  reduce  the  governing  partial  differential  equa- 
tion into  an  approximately  equivalent  ordinary  differential 
equation  will  be  sought.  Extension  of  the  proposed  stability 
theory  to  the  ignition  problem  will  be  continued;  attempts 
made  so  far  were  inconclusive.  Extension  of  both,  stability 
theory  and  computer  simulated  tests,  in  the  high  pressure  ran- 
ge (above  40  atm)  and  in  the  low  pressure  range  (below  5 atm) 
will  be  done.  Very  likely,  the  subatmospheric  pressure  range 
will  require  a careful  choice  of  the  flame  model,  including 
its  steady  state  aspects.  Studies  on  the  self  - sustained 
oscillating  regime  will  be  pursued. 

Experimentally,  data  on  combustion  vs  extinction  will 
be  furnished  both  from  a specially  designed  strand  burner  and 
a properly  modified  shock  tube.  This  should  allow  at  least 
a qualitative  comparison  between  predicted  and  experimentally 
observed  dynamic  stability  boundaries.  Better  data  are  expected 
from  the  laser  doppler  velocimetry  apparatus. 


